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&I Abstract 



< 



This study deals with asymptotic models for the propagation of one-dimensional internal 
w~) waves at the interface between two layers of immiscible fluids of different densities, under the 

,— ( rigid lid assumption and with a flat bottom. We present a new Green-Naghdi type model in 

the Camassa-Holm (or medium amplitude) regime. This model is fully justified, in the sense 
^_ ' that it is consistent, well-posed, and that its solutions remain close to exact solutions of the full 

'~j Euler system with corresponding initial data. Moreover, our system allows to fully justify any 

■^U( well-posed and consistent lower order model; and in particular the so-called Constantin-Lannes 

(-H approximation, which extends the classical Korteweg-de Vries equation in the Camassa-Holm 

■^-J regime. 



s 



1 Introduction 



'^ 1.1 Presentation of the problem 

l^L In the present paper, we study the propagation of internal waves in a two-fluid system, which 

iy-> consists in two layers of immiscible fluids of different densities, under the only influence of gravity. 

^^ The domain of the two layers is infinite in the horizontal space variable (assumed to be of dimension 

_J d = 1) and delimited above by a flat, rigid lid, a below by a flat bottom. Moreover, we assume that 

f~^ the fluids are homogeneous, ideal, incompressible and irrotational. We let the reader refer to [22| , 

C^ and references therein, for a good overview of the ins and outs concerning density-stratified fluids 

^~~^ in oceanography, and the relevance of our setup as a model for such system. 

^ The governing equations describing the evolution of the flow under the aforementioned configu- 

; ""j ration may be reduced to a system of two evolution equations located at the interface between the 

rS two layers (following a strategy initiated in the water-wave configuration in [15p9| , and achieved in 

C^ the bi-fluidic case in [8 ), named full Euler system. However, the study of this system is extremely 

challenging. In particular, the well-posedness of the Cauchy problem has been answered satisfac- 
torily (that is, with an existence of solutions on a time scale consistent with physical observations) 
only recently; see ^28] . 

Under these circumstances, a great deal of interests has been drawn to asymptotic m.odels, in 
order to predict accurately the main behavior of the system, provided some parameters describing 
the domain and nature of the fiow are small. Parameters of interests include 

dl a di pi 

A^ rfl d2 P2 



*IRMAR - UMR6625, Univ. Rennes 1, CNRS, Campus de Beaulieu, F-35042 Rennes cedex, France. | vincent.| 
duchene@univ-rennesl.fr 

TMathcmatiqucs, Faculte des sciences I et Ecole doctorale des sciences et technologic, Universite Libanaise, 
Beyrouth, Liban and CRAMS: Center for Research in Applied Mathematics and Statistics, AUL, Lebanon. 
|sJsrawi83@hotmaitcom 

fMafEematiques, Faculte des sciences I et Ecole doctorale des sciences et technologic, Universite Libanaise, Bey- 



routh, Liban. [r^lhoukSuLedu.lb, 



A new Green-Naghdi model in the Camassa-Holm regime April 17, 2013 



where a is the maximal vertical deformation of the interface with respect to its rest position; A 
is a characteristic horizontal length; di (resp. ^2) is the depth of the upper (resp. lower) layer; 
and pi (resp. P2) is the density of the upper (resp. lower) layer. Mathematically speaking, fj, 
and e measure respectively the amount of dispersion and nonlinearity which will contribute to the 
evolution of internal waves. It would be quite tedious to make an attempt at acknowledging every 
work on this aspect, but let us introduce some earlier results directly related to the present paper. 

Shallow water (/i <C 1) asymptotic models for uni-dimensional internal waves have been de- 



rived and studied in the pioneer works of 30 31133 . More recently, Choi and Camassa developed 



models with weakly (e = 0{fi)) and strongly (e ~ 1) nonlinear terms, respectively in 11 12 with 



horizontal dimension d — 2. They obtain bi-fluidic extensions of the classical shallow water (or 
Saint- Venant [I6]), Boussinesq [9 10 and Green-Naghdi [20[|36] models. Similar systems have been 



derived in [34| (with the additional assumption of 7 « 1) and in |14| (using a different approach, i.e. 
making use of the Hamiltonian structure of the full Euler equations). However, the aforementioned 
results are limited to the formal level. Let us mention now the work of Bona, Lannes and Saut [8] 
who, following a strategy initiated in [5|6] in the water-wave setting (one layer of fluid, with free sur- 
face), derived a large class of models for different regimes, under the rigid-lid assumption and with 
flat bottom (see also l4] where a topography and surface tension is added to the system, and fvf^ 
where the rigid- lid assumption is removed). The models derived in these papers are systematically 
justified by a consistency result: roughly speaking, sufficiently smooth solutions of the full Euler 
system satisfy the equations of the asymptotic model, up to a small remainder. 

Yet the consistency is only one of the properties that an asymptotic model shall satisfy, for its 
validity to be ascertained. Indeed, it leaves two important questions unanswered: for a large class 
of initial data (typically bounded in suitable Sobolev spaces) 

1. (well-posedness) do the full Euler system, as well as the asymptotic model, define a unique 
solution on a relevant time scale? 

2. (convergence) is the difference between these two solutions small over the relevant time scale? 

Two different strategies may be achieved in order to answer these questions. In the water-wave 
setting, Alvarez-Samaniego and Lannes p] proved the well-posedness of the full Euler system, as 
well as its stability with respect to perturbations of the equation. As a consequence, any well-posed 
and consistent asymptotic model is fully justified, in the sense aforedescribed. This holds true in 
particular for the (quasilinear, hyperbolic) Saint- Venant model, Boussinesq models (whose well- 
posedness has been extensively treated in [35J) and Green-Naghdi models (whose well-posedness 
have been proved in [3 23J29]); see [26] Appendix C] for a reader's digest of the numerous known 



results concerning asymptotic models in the water-wave situation). Unfortunately, proving the 
stability of the bi-fluidic full Euler system is yet out of reach, so that an other approach has to be 
used. In that case, one needs to obtain the stability of the asymptotic model itself, in order to deduce 
its full justification from a consistency result. Such a property holds true for the bi-fiuidic shallow- 



water system (see 21 ), and has been obtained in 18 for a class of Boussinesq-type systems. To our 



knowledge, no similar result is known concerning bi-fluidic Green-Naghdi models. Recall that Green- 
Naghdi models consist in higher order extensions of the shallow water equation, thus are consistent 
with precision 0{p.'^) instead of ©(/i), and allow strong nonlinearities (whereas Boussinesq models 
are limited to the long wave regime: e = 0{p)). Consequently, they contain nonlinear dispersive 
terms, which complicates their study, and in particular energy estimates. Finally, let us mention 
the recent work of Xu [38| , which studies and rigorously justify the so-called intermediate long wave 
system, obtained in a regime similar to ours: e ^ y^, but S ^ ^Ji. 

In this work, we present a new Green- Naghdi-type model in the Camassa-Holm (or medium 
amplitude) regime, e = 0{yj~p), for the propagation of internal waves. More precisely, the regime 
of validity of our model is, with fixed < /^max, M, Snun, Smax < 00: 

< A* < /"max, < e < Uim{My/JI, 1), S^in < S < ^niax, < 7 < 1 . 

Our model is fully justified: we prove that the full Euler system is consistent with our model, and 
that our system is well-posed (in the sense of Hadamard) in Sobolev spaces, and stable with respect 
to perturbations of the equations (so that the convergence property follows) . 
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Let us emphasize that in addition to its relevance as an asymptotic model, our system offers 
an important tool for the justification of other models. Indeed, it suffices to check that a given 
approximate solution solves our system up to a small remainder, to ensure that it is truly close to 
the solution of our Green-Naghdi model, and therefore to the solution of the full Euler system. 

Such strategy has been used in particular in [T] in order to rigorously justify the historical 
Korteweg-de Vries equation as a model for the propagation of surface wave in the long wave regime 
(with a Boussinesq model as the intermediary system); and this result has been extended to the 



bi-fluidic case in 18 . Higher order models in the Camassa-Holm regime have been introduced and 



justified in the sense of consistency in 13 in the water-wave case, and in 19 in the bi-fluidic case 



We are therefore able to fully justify the latter model, in the sense described above. 

1.2 Organization of the paper 

The paper is organized as follows. 

Section 2: The full Euler system. 

Section 3: Main results. 

Section J^: Construction of our model. 

Section 5: Preliminary results. 

Section 6: Linear analysis. 

Section 1: Proof of existence, stability and convergence. 

Section 8: Full justification of asymptotic models. 

Appendix A: Product and commutator estimates in Soholev spaces. 

We start by introducing in Section 2 the non-dimcnsionalized full Euler system, describing the 

evolution of the two-fluid system we consider. 

In Section 3, we present our new model, and we announce the main results of this paper. 

This asymptotic model is precisely derived and motivated in Section 4. 

Sections 5 and 6 contain the necessary ingredients for the proof of our results, which arc completed 

in Section 7. 

In Section 8, we explain how our system allows us to justify any well-posed and consistent lower 

order model. 

Finally the title of Appendix A is self-explanatory. 

We conclude this section with an inventory of the notations used throughout the present paper. 

Notations In the following, Co denotes any nonnegative constant whose exact expression is of no 

importance. The notation a < 6 means that a < Cq b. 

We denote by C(Ai, A2, . . . ) a nonnegative constant depending on the parameters Ai, A2,. . . and 

whose dependence on the Xj is always assumed to be nondecreasing. 

We use the condensed notation 

A, = Ss + (a),>,, 

to express that As ^ Bg li s < s and As — Bg + Cg if s > s. 

Let p be any constant with 1 < p < 00 and denote LP — LP{R) the space of all Lebesgue- measurable 

functions / with the standard norm 



\f\L. = {jjfixWdxy^'<^. 

The real inner product of any functions /i and /2 in the Hilbert space L^ (M) is denoted by 



(/l,/2)= / flix)f2{x)dx. 

Jm. 
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The space L°° = L°° (M) consists of all essentially bounded, Lebesgue-measurable fmictions / with 
the norm 

I/Iloc = esssup|/(a;)| < 00. 

We denote by W^'°° = W^'°°{R) = {/, dxf G L°°} endowed with its canonical norm. 
For any real constant s > 0, H'^ = H'^{R) denotes the Sobolev space of all tempered distributions 
/ with the norm 1/1//= = |A''/|i2 < 00, where A is the pseudo-differential operator A = (1 — 9^)^/^. 
For a given /i > 0, we denote by Hf^~^^{M.) the space _ff'*+^(IR) endowed with the norm 

I |2 _ I |2 I P 

I ■ \h^+^ = I ' Ih= + n ' Ih=+i • 

For any function u = u{t, x) and v{t, x) defined on [0, T)xM. with T > 0, we denote the inner product, 

the LP-norm and especially the L^-norm, as well as the Sobolev norm, with respect to the spatial 

variable X, by {u,v) = (u(t, •), w(i, •)), \u\lp = \u{t,-)\LP, \u\l2 = |u(i, •)|i2 , and \u\h- = |u(i, •)|h=, 

respectively. 

We denote L°°{[0, T); H^{M.)) the space of functions such that u{t, •) is controlled in H'^ , uniformly 

forie [o,r): 



\u\ 



L°°([0,T);ff=(R)) 



eSSSUp |u(t, •)|/f-a < 00. 

te[o,T) 



Finally, C'^(M) denote the space of fc-timcs continuously differentiablc functions. 
For any closed operator T defined on a Banach space X of functions, the commutator [T, /] is 
defined by [T, f]g = T{fg) — fT{g) with /, g and fg belonging to the domain of T. The same 
notation is used for / an operator mapping the domain of T into itself. 



2 The full Euler system 

We recall that the system we study consists in two layers of immiscible, homogeneous, ideal, in- 
compressible fluids only under the influence of gravity (see Figure nj. We restrict ourselves to 
the two-dimensional case, i.e. the horizontal dimension d = 1. The derivation of the governing 
equations of such a system is not new. We briefly recall it below, and refer to IlllslflQ] for more 
details. 

Z 



. 
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Figure 1: Sketch of the domain and governing equations 

We assume that the interface is given as the graph of a function C(t, a;) which expresses the 
deviation from its rest position (a;,0) at the spatial coordinate x and at time t. The bottom and 
surface are assumed to be rigid and flat. Therefore, at each time i > 0, the domains of the upper 
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and lower fluid (denoted, respectively, fij ^^^ ri|), are given by 



n\ 



^2 



{(2:,z)eMxM, C{t,x) < z < di}, 
{ {x,z) eRxR, -d2 < z < C{t,x)}. 

We assume that the two domains are strictly connected, that is there exists h > such that 

di-({t,x)>h>0, and d2 + C{t,x) > h > 0. 

We denote by (pi, vi) and (p2, V2) the mass density and velocity fields of, respectively, the upper 
and the lower fluid. The two fluids are assumed to be homogeneous and incompressible, so that the 
mass densities pi, p2 are constant, and the velocity fields Vi, V2 are divergence free. As we assume 
the flows to be irrotational, one can express the velocity field as gradients of a potential: v^ — V0i 
{i = 1,2), and the velocity potentials satisfy Laplace's equation 



di 



di<j>^ = 0. 



The fluids being ideal, they satisfy the Euler equations. Integrating the momentum equations 
yields the Bernoulli equation, written in terms of the velocity potentials: 



dt(t)i 



2"^- 



P 

P 



--gz innl (i = l,2). 



where P denotes the pressure inside the fluids. 

From the assumption that no fluid particle crosses the surface, the bottom or the interface, one 
deduces kinematic boundary conditions, and the set of equations is closed by the continuity of the 
stress tensor at the interface, which reads 



lP{t,x)l = lim(p(i,2;,C(i,x)+e) - F(t, a;, C(i, a;) - e) ) 
■7^ 



Tk{at,x)), 



where k(C) — — ^a-l —r= ^ dxC] denotes the mean curvature of the interface and a the surface 

tension coefficient. 

Altogether, the governing equations of our problem are given by the following 



(2.1) 



where n denotes the unit upward normal vector at the interface. 

The next step consists in nondimensionalizing the system. Thanks to an appropriate scaling, the 



dl(j), + dl(j), = 
8,(1)1 = 


~ gz 


in fi*, i = 1,2, 
in VL\, i = 1,2, 
on Ft = {(x, z), z = di}. 


dtC = Vl + |5xCP5„<^i 
9,02 = 
lP{t,x)l =. ak{Q 


= ^l + \d^C?dnh 


onT = {{x,z),z = Q{t,x)} 
on Fb = {{x,z),z = -^2}, 
on F, 



two-layer full Euler system (2.1) can be written in dimensionless form. The study of the linearized 



system (see 28 for example), which can be solved explicitly, leads to a well-adapted rescaling. 



Let a be the maximum amplitude of the deformation of the interface. We denote by A a 
characteristic horizontal length, say the wavelength of the interface. Then the typical velocity of 
small propagating internal waves (or wave celerity) is given by 



Co 



(P2 - Pi)did2 



V ^2^1 + ^1^2 

Consequently, we introduce the dimensionless variablefjj 



z 
^1 



X 

A' 



t = 



Co 
A 



t. 



-"^We choose d\ as the reference vertical length. This choice is harmless as we assume in the following that the two 
layers of fluid have comparable depth: the depth ratio 5 do not approach zero or infinity. 
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the dimcnsionlcss unknowns 

a 
as well as five dimensionless parameters 



Pi 
P2' 



4)i{i,x,z) 



dl 



aXcQ 



dl 



i{t,x,z) (i = l,2), 



dx- 



A2' 



Bo 



g{p2 -pi)A^ 



We conclude by remarking that the system can be reduced into two evolution equations coupling 



Zakharov's canonical variables 15 , 39 , namely (withdrawing the tildes for the sake of readability) 



the deformation of the free interface from its rest position, ^, and the trace of the dimensionless 
upper potential at the interface, ^, defined as follows: 

ip = c/)i{t,x,C{t,x)). 

Indeed, 0i and 02 are uniquely deduced from (C, ip) as solutions of the following Laplace's problems: 



(2.2) 



(2.3) 



{^idl + d^) (J3i=0 in Oi = {{x,z) e R^, eC{x) < z < 1}, 



9.01 = 

01 =^ 

( t^di + dl ; 

dn(t>2 = 9„(/)i 

9.02 = 



on Tt = {{x,z) e 



= 1}, 



onr = {{x,z) effi^ z = tC}, 

i2 = inf72 = {(x,z) eM^^ -i <z<eC}, 



onT, 

onrb = {(.T,z) eM^ 



More precisely, we define the so-called Dirichlet-Neumann operators. 



Definition 2.1 (Dirichlet-Neumann operators). Let ( € H*°^^{M.), ig > l/^i such that there exists 
h>0 with hi = I- eC,>h>Q andh2 = \ + eC, > h > {), and let ij) e Lf^^{M.),dxi' e H'^/'^{R). 
Then we define 



G''^ =G''[e(]yj = ^1 + Mk9,CP(9„0i) !.=< = -/ie(a.C)(9.0i) U=.c + (5-<^i) U=< ^ 

^M^ =i/A<'*[<]^ = d^{cj)2U=ec) = (a,02)|.=<+e(9,C)(5.02)|.=eC, 

where 0i and 02 are uniquely defined (up to a constant for 02 J as the solutions in _ff^(M) of the 



Laplace 's problems { 2.2 ) -( 2.3 1 



The existence and uniqueness of a solution to (2.2)-(2.3), and therefore the well-posedness of 



the Dirichlet-Neumann operators follow from classical arguments detailed, for example, in 126]. 
Using the above definition, and after straightforward computations, one can rewrite the nondi- 



mensionalized version of (2.1 ) as a simple system of two coupled evolution equations, namely 
1 



(2.4) 



dtC - -G^^P = 0, 



a,(i?''^V-79.^) + il + S)d,C + |5,(|i/^'^^|2-7|9,^|' 



fied^M'^'^ 



-f + 5 d^{k{e^C)) 



Bo 



ey^ 



where we denote 



N^'^ = 



2{\ + ^i\e^M'') 



We will refer to (2.4) as the full Euler system, and solutions of this system will be exact solutions 
of our problem. 



April 17, 2013 



Vincent Duchcnc, Samcr Israwi, Raafat Talhouk 



3 Main results 

We now present our new Green-Naghdi model, that we fully justify as an asymptotic model for the 



full Euler system (2.4), in the so-called Camassa-Holm regime, that is when e = 0{y/Ji),n <^ 1. 
More precisely, we assume that there exists < M, /Zmax, ^min, ^max < oo such that the dimensionless 
parameters (e,/i, (5, 7) belong to the following set: 
(3.1) 

VCH = {(/i,e,(5,7) : < ^ < /imax, < e < min(7\f ^, 1), S <E ((5min, ^max), < 7 < 1 }. 

All the estimates are now understood uniformly with respect to (/x, e, S, 7) £ Vch- 
Our new system is: 



(3.2) 



T[eC] {dtv + e^vd^v) + (7 + 6)qM)dxC 



where hi = 1 — e(^ (resp. /i2 = | + e^) denotes the depth of the upper (resp. lower) fluid, and v is 
the shear mean we/ocJ<j^ defined by 

(3.3) -G-^K]^ = -Mir^')- 

The operator 'X is as follows: 

(3.4) 1[e(:]V = qM)V - tii^d^[q2{eOd.v), 



with qi{X) = 1 + KiX (i = 1,2) and ki,K2j'^ are constants displayed in (4.9) and (4.121, later on. 
This model is fully justified by the following results: 

Theorem 3.1 (Consistency). Let U^ = {(^jtp^) be a family of solutions of the full Euler sys- 
tem (2.4) such that there exists Co,T > with 



C" 



lL°°([0,T);ff=+2) 



\dtC'\ 



J A-\\f) ,/,P 



L-([0,T);ff»+V) + ll^*^-^iL-([0,r);ff=+l) - <^0' 



for any s > So + 1/2, sq > 1/2, and uniformly with respect to (/i, e, ^, 7) = p G Vch; see (3.1) 
Moreover, assume 



(HI) 



3hoi > such that hi = 1 - eC*" > ^loi > 0, ^-2 = 7 + eC > /loi > 



Define v'^ as in (3.3). Then (C,^) satisfies (3.2|, up to a remainder R, hounded by 

II^IIl~([o,t);..=) < (M^+Bo-i)C, 
with C — C(/iq]^ , M, /irnax, 5^5 (^max, Cq) , Uniformly with respect to the parameters p G Vch ■ 

Our system is well-posed (in the sense of Hadamard) in the energy space X'^ = H''{M.) x iJ''+^(M), 
endowed with the norm 

provided the following ellipticity condition (for the operator T) holds: 

(H2) 3ho2 > such that inf (1 -I- eK2C) > /102 > ; inf (1 + ekiC) > /io2 > 0. 

xeR xeR 

^v is equivalently defined as v = U2 — 7M1 where ui, U2 are the mean velocities integrated across the vertical layer 
in each fluid: Si (t,x) = ^^(f^^.) f^c^t.x) dx<f>i{t,x,z) dz and U2{t,x) = ^^(f^j;) f^i'^'dx<f>2{t,x,z) dz. 
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Theorem 3.2 (Existence and uniqueness). Let sq > 1/2, s > sq + 1, and let Uq — (Co,i'o)^ G X'' 
satisfy (|H1|),(H2 1. Then there exists a maximal time Tmax > 0, uniformly hounded from below 



with respect to p G VcH , such that the system of equations ( |3.2[ ) admits a unique solution U = 
(C,«)^ e CO([0,T„iax/e);^'^)nCi([0,T„iax/e);^""') with the initial value {C,v) |t=o = (Co,wo) and 
preserving the conditions ( |Hl| , pI2| (with different lower bounds) for any t € \Q,T^g_^/e). 

Moreover, for any < T < T^ax, there exists Cq^Xt = C(/ioi\ /io2\^max,^~?n, Af,r, |t/o|^J, 
independent of p £ TcH, such that one has the energy estimate 



W 0<t< 



T 



|c/(t,-)L, + |a,c/(i,-)L._i <Coe 



eXTt 



//^max < oo, one has 



\U{t,-)\x.^ 



oo as t 



T„ 



or one of the two conditions ( |Hl[ ), pr[2| ceases to be true as t — > Tn^a^/e. 

Theorem 3.3 (Stability). Let s > sq + I with sq > 1/2, a nd le t [/o,i = (Co.i; ^o,i)^ G X'^ and Uq,2 
(Co. 21^0,2)^ G X'^^^ . Under the assumptions of Theorem 3.2, let Uj be the solution of system (3.2) 



with Uj \t=o = Un.j-Then there exists T,X,Cq = C{hgl,hg2 , <^max, ^,„in' ^^' |C^o,i|jf., |t^o,2U=^+i) > 
such that 

Vt€[0,^], \iUi-U2){t,-)\xs <Coe'^^'\Ui,o-U2,o\^^. 

Finally, the following "convergence result" states that the solutions of our system approach the 
solutions of the full Euler system, with as good a precision as ^ (and Bo^^) is small. 

Theorem 3.4 (Convergence). Let p e Vch and U^ = (C'',V''°)^ G iJ'*+^, N sufficiently large, 
satisfy the hypotheses of Theorem 5 in l^28r\ as well as (H1),(H2). Then there exists C,T > 0, 
independent ofp, such that 

• There exists a unique solution U = (CV')^ lo the full Euler system (2.4), defined on [0,T] 



and with initial data (C , "0 ) (provided by Theorem 5 in '28^); 



• There exists a unique solution Ua = {Ca^^a)'^ to ou r new model (3.2 I, defined on [0,T] and 
with initial data {C^,v^)'^ (provided by Theorem 3.2); 



With V = v[C,,ij], defined as in (3. 3), one has 



|(C,«) - (Ca,«a)L^([o,j,].^,) < Cifl' +B0~')t. 

Remark 3.5. The above proposition is valid for time interval t G [0, T/e] with T bounded by below, 
independently 0/ p G VcH, provided that a stronger criterion is satisfied by the initial data ; see 
criterion (5.5) and Theorem 6 in ^281. 

Our new model allows to fully justify any well-posed asymptotic model, consistent with our 



model (3.21 



Proposition 3.6. Consider (S) a system of equations such that 

• The Cauchy problem for (S) is well-posed in X^ , r sufficiently large. 

• For U^ = (C",!;*')^ G H^^^ , N sufficiently large, the solutions of (S) satisfy our model (3.2), 



up to a remainder R of size 0{i) in L°°{[0,T]; H" x H'^). 
Then under the assumptions of Theorem\3.4\ the difference between the solution of the full Euler 



system (2.4), U = (CV")^; o,''^d the solution of the asymptotic model (S) with corresponding initial 
data, Ua = (Ca,Wa) , is estimated as follows: 

I (C, V[Q, 0]) - (Ca, «a) L^([o_y].^.) < C{i + /? + B0-l)i. 

This procedure of full justification is precisely described in section|8J and applied to the so-called 
Constantin-Lannes decoupled approximation model. 



^in particular, it satisfies a stability criterion, which in the shallow wat er c onfiguration (/^ <S 1), can be roughly 



expressed as T ; 



^^^2° j^ YTL \o is sufficiently small ; see section 5.1.3 of pSi. 
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4 Construction of the model 

This section is dedicated to the construction of the model we study. The key ingredient for construct- 
ing shallow water asymptotic models lies in the expansion of the Dirichlet-to-Neumann operators, 
with respect to the shallowness parameter, /i; see Proposition |4.1[ below. Thanks to such an ex- 
pansion, one is able to obtain the so-called Green-Naghdi model (for internal waves), displayed 



in (4.6). This model has been introduced by one of the author in [19|, and justified by a consis- 



tency result (Proposition 4.2): roughly speaking, any solution of the full Euler system satisfies the 
Green-Naghdi asymptotic model up to a small remainder, of size C(/x^ -|- Bo"^). n 

In a second step, we use the additional assumption of the Camassa-Holm regime, e = 0[^JJi). 
Simplified models, with the same order of precision may then be deduced. We use several trans- 



formations which allow to obtain a well-prepared model: system (4.13), presented in the previous 



section (3.2). The justification of this model, in the sense of consistency, is stated in Theorem 4.4 
The stronger results (well-posedness, stability, convergence) described in Section ^ are proved in 
subsequent sections. 

4.1 The Green-Naghdi model 

The following Proposition is given in 19 , extending the result of [8]. 



Proposition 4.1 (Expansion of the Dirichlet-Neumann operators). Let s > sq + 1/2, sq > 1/2. 
Let i) he such that d^ij) G if'*+"/2(M), and C, £ H''+'^/'^{M.). Let /ii = 1 - eC and /i2 = 1/5 + eC, such 
that there exists h > with hi,h2 > h > 0. Then 

(4.1) |-G^^ - d,{h,d,i;)\^^ < fiC,, 

(4.2) |-G^V - d,{h,d,ij) - fi\dl{hldli;)\^^ < ^2^3, 

l-l O 

(4.3) |i?^''^+ 5^9,^1^, < fiCo, 



(4.4) lif'-V + j^a,^ - Jl-d.,[hld4^d,^) - hldl^p) 1^., < ^^^ c^. 



with Cj — C{j^, /imax, y^i ^max, | C| j:/s+3/2+3 j l^xi^] fjs+5/2+j) ■ The estimates are uniform with respect 
to the parameters e € [0, 1], ^ G [0, ^max], S G ((5min, ^max)- 



Plugging these expansions into the full Euler system ( 2.4 ), and withdrawing 0{iJ, ) terms, imme- 
diately yields an asymptotic Green-Naghdi model. This model is presented in [19| and justified in 
the sense of consistency. However, such a Green-Naghdi model is only one of the variety of models 
which satisfy such a property. In the following, we decide to introduce an equivalent model, using 
as unknown (instead of ^, the trace of the lower velocity potential at the interface) the shear mean 
velocity, defined by 

(4.5) V = U2 — 7Mi, 

where Ui, U2 are the mean velocities integrated across the vertical layer in each fiuid: 

ui(t,x) = T-— — r / dx(t>i{t,x,z) dz, and U2{t,x) = -—t- — r / dx(l>2it,x,z) dz. 

hi[t,x) J^^(^t,x) h2it,x)J_i 



*We neglect the surface tension term in our asymptotic models, that is we set Bo~ = 0. As a matter of fact, the 
surface tension effects are necessary for the full Euler system l |2.4[ l to be well-posed. However, the surface tension 
is very small in practice, and does not play any role in the behavior of our asymptotic models. The idea behind 
this apparent paradox is that the surface tension helps to regularize the high frequencies of the full Euler system, 
whereas the solutions of our models are located in low frequencies, where surface tension effects are negligible. We 
refer to [28] for an in-depth analysis of this phenomenon. 
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Equivalently (as shown in fTo"), one has 

1 






Such a choice has been used in 11 12 for example, and present at least two benefits. First, the 



equation describing the evolution of the deformation of the interface is an exact equation, and not 
a 0{iJ?) approximation. Indeed, one has immediately from the full Euler system (2.4): 



dtC 






hih2 



7/12 



What is more, the system obtained using mean velocities have a nicer behavior with respect to its 
linear well-posedness, thus one can expect nonlinear well-posedness only for the latter. 

Altogether, several technical but straightforward computations yield the following Green-Naghdi 
model (without surface tension, Bo""'^ = 0) 
(4.6) 



dtiv + ^iQ[h,,h2]v\ + {-/ + 6)d.,C + ^^-((^^^7^2 1^1") = ^^e^,{n[h,,h2]v) 



where we define: 



Qlh.Mv - ^.M'^Mj^j) + ih2d.{hlo4-l 



h2 V 



-7/12 



hi V 



n[hiM]v = ^( [h2d:c{ J"\ 1, ))' - i{hid.,{ " 



'hi +7/12 



/l2 V 



hi + 7/12 



1 V 



hi 



^dJhld4- 
h2 \ ^hi 



hi V 



1- 



e<,.(.3<,.(^ 



/i2 y 



3/11+7/12 l/i2 V ^/ii+7/i2 / hi \ ^hi+jh2- 

This model has been derived in |19 and justified in the sense of consistency, as follows. 



Proposition 4.2. Let U^ = (('P,'!/''^) be a family of solutions of the full Euler system (2.4 1 such 
that such that there exists Co,T > with 



e 



\dtC\ 



-i. + P.i^'l 



\dtd,^^\ 



.+ 9 < Co, 



lL°°([0,r);ff°+3) II ^ llL°°([0,r);ff°+2) II ^^ H L~ ([0,T);H'' + ^ ) ' W"'"-^'^ ^^ L°°{[0,T);H''^i) 

for any s > sq + 1/2, sq > 1/2, and uniformly with respect to (/i, e, 5, 7) = p G Vsw ■ 
'Psw = {(a«, e,^,7) : < ^i < /x^ax, < e < 1, S^^in < S < S,^^^, < 7 < 1 }. 
Moreover, assume that there exists h > such that 

/ii = 1 - eC" > /i > 0, h2 = l+eC >h>0. 





Define v'^ as in (4.5) or, equivalently, by 



fi \hi + 7/12 



Then {C,v) satisfies (4.6), up to a remainder R, bounded by 



R 



L°°([0,T);ffO 



H.) < (/i'+Bo-i)C, 



with G = G{j^, Hynax, -f^ , 5max, Gq) , Uniformly with respect to the parameters p G Vsw- 
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Remark 4.3. In ^191, the author works with Bo = 0. However, it is straightforward to see that 
the surface tension term is estimated as follows: 

j + Sd^kje^O) 
II Bo e^ llL~([o,T);ff=) - ^o ^' 

with C = C((5max, C rooffn TV^s+^O' ''uniformly with respect to the parameters p G Vsw- 

4.2 Our new model 

As announced in the introduction, the present work is Hmited to the so-called Camassa-Holm 

regime, that is using additional assumption e = 0{y/Ji). More precisely, we assume that there 

exists < M, /imax, (^min, (^max < oo such that the dimensionless parameters (^, e, 5, 7) belong to the 

following set: 

(4.7) 

VcH = {(^i, e,(5,7) : < ^ < ^,„ax, < e < min(A/^, 1), Jmin < ^ < (^max, < 7 < 1 }. 

All the estimates are now understood uniformly with respect to (/i, e, 5, 7) e Vch- 
In particular, one can easily check that the following approximations hold: 

Q[huh2]V = -lydlv - e^^ {ip - a)vd'',C + [a + 2l3)d,iCd^v) - ^Cdlv) + 0{y), 



TZ[h,,h2]V = a(^{d,vf + ^vdlv] + 0{^) 



with 

(4-8) V = — - — '—- , a = 7 ^ and ^ = ^ ^ ^,, 



Plugging these expansions into system (4.6) yields a simplified model, with the same order of 
precision of the original model (that is 0{iJr)) in the Camassa-Holm regime. However, we will use 
several additional transformations, in order to produce an equivalent model (again, in the sense of 
consistency), which possess a structure similar to symmetrizable quasilinear systems, which allows 
the study of the subsequent sections. 

The first step is to introduce the following symmetric operator, 

1[e(:]V = qM)V - fii^d,(^q2{eOd,v), 
where qi{eC) = 1 + '^jeC (* = li 2) and ki, K2 are constants to be determined, so as to write 
qi{eC)dt(v + ^J-Q[hi,h2]v] = 1[eQdtV + higher order terms. 
More precisely, one can check 
1[eQdtV - qi{eC)dt(v + fJ,Q[hi,h2]v) = - fieiyK2dxUdxdtv] + iiev Ki^dldtv 



+ ti-^qii 



{eC)^-dt[iP - a)vdlC + ia + 2p)dxiCdxv)-^Cdlv 



Now, we develop the latter term, making use of the fact that at first order, system (4.6 1 yields (in 
the Camassa-Holm regime (|4.7|)) 



dtv ^ -{^ + 6)d,C + 0{^) and dtC ^ ^d^v + 0{^). 

7 + (5 
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It follows 



+ fMhi,h2]v^ = nei^{-f + d)(^K2d^(^Cdlc) - '^iC^^c) 



^(7 + -5) 



^'3^ ((/3 - c^)d.CdlC + (a + 2/3)9. (Ca^C) - PCd^c) 



He 



^(^{13 - a)vdlv + ia + 2l3)dAid,v)^) - P{d^v)idlv)) + 0(/i'). 



It becomes clear, now, that one can adjust ki^ K2 £Rso that all terms involving <^ and its derivatives 
are withdrawn. More specifically, we set 



(4.9) 



Ki = 



-f + 5,^^ _N o'^^^T fl^T .. {l + S)[3 ^S^~l 



-{2/3 -a) = 2— -^-<5-— -, K2 = 
Siy 7 + 1 + 7c) 



= 3 



7 + ^' 



and one obtains consequently 
(4.10) 



However, one of the remaining terms in (4.10), as well as in dx\TV[hi^h2]vj, involves three 
derivatives on v. In order to deal with these terms, we introduce T[eC](e<?w9a;'y) where, again, <; G M 
is to be determined. More precisely, one has 

%[eC\{ei,vdxv) + ^ieqi{eC,)dx{Tl[hi,h2]vj 

= ei,qi{eC,)vdxV - ^iel'(;^x\Q2{eC)^x{v^xv)]+^^€ql{€C)a^x(-{^xvf + -v^lv 



In the Camassa-Holm regime (4.7), this yields 
(4.11) 



1[eC]iec;vdxv)+fieqi{eC)dx(n[hi,h2]v) ^ e^qM)vdxV+fiedx(^{^-'y';){dxvf + {'^-,^<;)vd^v)+0{fi^). 



(4.12) 



Combining (|4.10| with (|4.11|), one can check that if we set 

1-7 



2a -/3 



= 26- 



3iy "(7 + (5)(l+7(5) (j + S)^' 

then the following approximation holds (withdrawing 0(/i^) terms): 

1[eC]idtv + e(;vdxv) - qi{eC)dt(v + ^J.Q[hi,h2]v) + fieqi{eC)dx(TZ[hi,h2]v) 



P-a 



vdiv 



= e<;qi{eC,)vdxV - ne—dx{idxvf). 



When plugging this estimate in (4.6), and after multiplying the second equation by <?i(eC), one 
deduces that an equivalent model (that is, with the same precision for consistency) is 



(4.13) 



dtC + dx 



hih2 



v] ^ 0, 



hi +7/12 

T[eC] [dtv + e^vdxv) + (7 + 5)qM)d.C 
+ i9iK)^-((7g 



1+7/1 



:ih^\y\2 _ A^\2 



-lieladx{{dxvf), 



where we recall that 
(4.14) 



T[eC]V^ = qM)V - iiydx(q2{eOd.v), 



with qi{X) = 1 + KiX (i = 1, 2) and ki, K2, ^ are defined by ( |49| ) and \A.12\ . 

This is the system we study. We state below its validity as an asymptotic model for the full 
Euler equations in the Camassa-Holm regime, in the sense of consistency. This is Theorem |3.1[ in 
section [3l 
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Theorem 4.4. Let [/P = (C'', V''') be a family of solutions of the full Euler system (2.4) such that 
such that there exists Cq,T > with 



e 



lL°°([0,T);ff=+2) 



\dtC 



IL°°([0,T);H=+2) 



\d,^' 



lL-([0,T);H»+V) + ll^*^-^'llL-([0,T);ff=+l) " ^°^ 



for any s > sq + 1/2, sq > 1/2. and uniformly with respect to (/i, 6,^,7) = p G Vch; see (4.7) 
Moreover, assume that there exists h > such that 



hi 



l-eC >h>0, /i2 



1 



eC >h>0. 



Define v''' as in (4.5) or, equivalently, by 






hih2 



hi + 7/12 



Then {C,v) satisfies (4.13), up to a remainder R, bounded by 



R 



L~([0,T);ff=) 



< (^2+Bo-i) C, 



with C = C(t^, M, /ifjiaxi x^i (^maxj C'o), Uniformly with respect to the parameters p G Vch- 

Proof. Let U = {(,, ip) satisfy the hypotheses of the Proposition (withdrawing the exphcit depen- 



dence with respect to parameters p for the sake of readability). As a consequence of Proposition 4.2 



we know that (Cf) satisfies (4.6), up to a remainder Rq, bounded by 



with C = C(i,M,M: 



max, 3;;^,'5max,C'o), Uniformly with respect to (/[z,e, ^,7) G Vch 3 Vsw- The 
proof now consists in checking that all terms neglected in the above calculations can be rigorously 
estimated in the same way. 

The formal expansions can easily be checked. When turning to control the remainders in H'^ 
norm, the main issue we face is to obtain tame estimates which involve Ur- L .^ , controlled by h~^ . 
Let us describe how such a result can be obtained, with a simple example. 



Lemma 4.5. Let f,C^ L°° p| iJ*, with s > and hi — 1~ e^, with hi(e() > h > for any x G 
Then one has 

\lj\^. < C{h'\e\C\^J{\f\^, + e\C\^.\f\,^). 



Proof of the Lemma. We will use the identity 



hi 1 - eC 



< 



1-eC 



/• 



Using Moser's tame estimates (Lemma A.l ), one has 



'/ii 



f\ < \f\ 
< \f\ 



< 



/|^=- 



a function F € C° 



such that 



1-eC 

^eClH' 



\f\H^^ 



< 



1-eC 



l^=-|/| 



The only non-trivial term to estimate is now j^r^ \ fjs ■ Using that hi — l~e( > h > 0, we introduce 



F{X) = ^-^ 



a 1 - X > h > 0, 

iil-X <o. 



The function F satisfies the hypotheses of Lemma [A. ![ and one has 



< 



'1-eC 
The Lemma is now straightforward. 



H. = InOln^ < Ci\eC\^^,h~')\eC\^,. 



n 
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Let us apply this Lemma to the rigorous estimate of specific expansions. 

Lemma 4.6. Let f,C,e £°°([0, T);H^), with s > sq > 1/2, such that hi = 1 - €(,> h> Q for any 

(t,x) € [0,T) X M. Then one has 

, 1 



/ii 



/-/|^. < eCih-\e\C\^,)\C\^,\f\^, 



Proof. Let us first remark that the formal expansions are straightforward: 
1 . 1 



1 + 0{X) 



and 



i + x + oix-" 



for 



1-X ' ' \-X 

The rigorous estimate is obtained thanks to Lemma |4.5[ applied to the identities 

'/-/ = ?/ 



1^1 <1. 



h\ h\ 

Indeed, one has for s > sq > 1/2, 



f e^e 

and ^/-(l + <)/ = ^/- 



l^/-/|.^- - \fj\.-^ ^ l<l.=-|^/|.^- < ^Cih-\e\CU\f\.^, 



where we used Lemma A.l Lemma 4.5 and continuous Sobolev embedding, since s > sq > 1/2. 
The second estimate, concerning I -j^ f — (1 + eC)/|f7si is obtained identically. D 

Of course, the same results would follow when replacing hi by h2, or hi + 7/12 (with straight- 



forward adjustments). The expansions used throughout section 4.2 are more involved technically, 
but follow the exact same method. 

For the sake of brevity, we do not develop each of these results, but rather provide the exact 
estimate that one obtains. In order to simplify the presentation, we do not precise the dependence 
with respect ■^, M, //max, y^j <^max as p = (/i, e, S, 7) £ Vch in the following constant, but we only 
focus on the regularity required for Cw, and the accuracy of the estimates. 

One has 

\\dt{Q[hi,h2]V)-[-vdldtv-e^-^dt ((/3 - a)vdl(: + (a + 2/3)9,(C9.«) - Kdlv) ] |L^([o^^).^.) 



L=°([0,T);_f/=+3)' ln*^llL°°([0,T);ff=+2)' Ir llL=([0,T);_f/=+3)' lrllL~([0,T);ff»+2) 



with C{s + 3) = C(||C 
and 

\\d.{n[hi,h2]V) - d,[a{\{d^vf + \vdlv)] W^^^^^^^y^^.^ < eC{s + 3). 

Recall that e < M^/Jl. Then, since (C, v) satisfies ( |4.6[ ), up to the remainder Rq, one has 

1 



dtv + {-f + S)dM 



L~([0,T);ff-) 



dtC- 



^ + ^-^-^IL-([o 



.T);_ff=) - V7i'^(''^ + 3) + ||^o||ioo([o,T);ff=)- 



It follows that (4.11 ) is valid up to a remainder Ri, bounded by 

ll-Rill < M^C(s + 3)+^3/2||^^ 



L~([0,T);ff=)- 



Finally, {C,v) satisfies (4.61, up to the remainder Rq + Ri, and 

\\Ro + ^i|| < fJ- C'oj 
where we use that 

The estimate on v follows directly from the identity dx[ f^^\ v\ = —-G^'*^^ = dtC- The esti- 
mate on dtv can be proved, for example, following flTJ Prop. 2.12]. This concludes the proof of 
Theorem lO □ 
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5 Preliminary results 



In this section, we study the operator 'J[eC]j defined in (4.14): 

(5.1) %[eC\V = (1 + eniQV ~ ^iiyd, ((1 + eKaO^.V^) . 

with v,K,i,K2 are constants. In our setting, v,ki,K2 depend on the parameters 7,5; but in what 
follows, we use only that i^ > is bounded by below and i/ + |ki| + |k2| is bounded by above by 
constants depending only on Syam and 5a\ax 



1 + 7^ 
3(5(7 + 5) 



max 



ae. 



\^il\ + \n2\<C{5zl,5. 



mill' "max;, 



when we restrict to (/x, e, (5, 7) G Vch, as defined in (4.7), and that we recall below: 

VCH = {{p,e,S,j) : < /i < ^^max, < e < niin(Af^, 1), Snun < S < (5max, < 7 < 1 }. 

When no confusion is possible, we write simply T = 'J[eC]. In the following, we seek sufficient 
conditions to ensure the strong ellipticity of the operator T, independently of (^, e,(5, 7) G Vch, 
which will yield to the well-posedness and continuity of the inverse 'Z~^. 

As a matter of fact, this condition, namely (|II2[) (and similarly the classical non-zero depth con- 



dition, (HI )) simply consists in assuming that the deformation of the interface is not too large. For 
fixed C G L°° , the restriction reduces to an estimate on emax|C|r=o, with e^ax = min(A/.y/i„iax, 1), 



and (II1)-(II2) hold uniformly with respect to (/i, e,(5, 7) S Vch', see Lemma 5.1 below. 
Let us shortly detail the argument. Recall the non-zero depth condition 



(HI) 



3 hoi > 0, such that min (inf hi, inf /i2) > ^01 



where hi = 1 — e( and /i2 = j + e<^ are the depth of, respectively, the upper and the lower layer of 
ffiiid. It is straightforward to check that, since for all (/i, e, (5, 7) G Vch, the following condition 



c 



< min(l, ) 



is sufficient to define /iqi > such that (HI I is valid, independently of (/x, e,(5, 7) G Vch- Briefly, 
since e < e^ax, infa^eR /ii > l~e|C|icx, > l-emax|C|ioo, inf^eR /12 > j - e|C|ic=o > j " emax|C|ioo- 



Note that conversely, for (HI) to be satisfied for any (e,^, 5, 7) G Vch, then one needs 

1 . 



c 



< maxfl, 



^min 



Indeed, one has for e = e,„ax, CmaxC = 1 - hi < 1 - /iqi and -e„iaxC = j - ^2 < j - ^oi- 
In the same way, we introduce the condition 

(H2) 3 ho2 > 0, such that inf (1 -I- eKzC) > ^02 > ; inf (1 + eKiC) > /io2 > 0. 

xeR xeR 

As above, such a condition is a consequence of a simple smallness assumption on e C roo- More 
precisely, one has the following result. 

Lemma 5.1. Let C G L°° and Cmax = min(Af ^/imax, 1) be such that there exists ho > with 

max(|K2|,|Kl|,l,(5,„ax)emax|C|roo < ^ ^ ^q < 1. 



Then there exists /ioi,^o2 > such that (H1|-(H2| hold for any (/i, e,(5, 7) G Vch 



In what follows, we will always assume that (HI I and (H2) are satisfied. It is a consequence of 



our work that such assumption may be imposed only on the initial data, and then is automatically 
satisfied over the relevant time scale. 
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Before asserting the strong eUipticity of the operator T, let us first recall the quantity | • |//i, 
which is defined as 

and is equivalent to the i/^(M)-norm but not uniformly with respect to ii. We define by Hl^{R) the 
space i?^(M) endowed with this norm. 



Lemma 5.2. Let ( e L°°(R.) be such that (H2) is satisfied. Then the operator 



TK] : H^iR) -^ iHl{R)r 
is uniformly continuous and coercive. More precisely, there exists Cq > such that 



(5.2) 
(5.3) 



(Tu,u) < co\u\hi\v\hI: 

{1u,u) >-\u\j,^ 
Co f 



with Co = C(/io2 ,1^,1' M«iUK2|,e|C|i^)- 
Moreover, the following estimates hold: 

(i) Letso>lands>O.If(e iJ^«(M) n iJ"(M) and u G H^+^{R) and v e H^{R), then: 

(5.4) \{A^1[eC]u,v)\ < Co((l + e|C|^.„)h|^.+i+(e|C|H.|"l//;o+i),>,J|Hffi > 
(it) Let So > 5 ands>0. If ( E H''«+^ D iJ*(E), u E H^iR) and v E H^{R), then: 

(5.5) |([A^T[eC]]u,«)| < e Co (|C|^.o+ih|^. + (|C|^.hlH;"+0.>so+i) I^Ih; ' 

where Cq = C(/lo2^ ^j ^^"^j ki|) I'^iD- 
Proof. Let us define the bilinear form 

a{u, w) = ( Tw , «) = ((! + eKi()u , w ) + ^/i( (1 + eK2C)dxU , d^v ) , 

where ( • , • ) denotes the L^-based inner product. It is straightforward to check that 

|a(u,w)| < sup |1 + eKiC|( u , w ) + v/isup |1 + eK2C|( dxU , d^v ), 

xeR xeR 



so that (5.2) is now straightforward, by Cauchy-Schwarz inequality. 

The iJ^(M)-coercivity of a(-, •), inequality (5.31, is a consequence of the condition (H2) 



a{u, u) — ( Tu , u ) = / (1 + eKif^)!?*! dx + v^i / (1 + eK2C)l''^2;| dx 
> /io2min(l,i/)|u|^i. 



Let us now prove the higher-order estimates of the Lemma, starting with the product estimates. 
One has 

( A'Tu , v)^{A'{{l + en.Ou} , v ) + iy^,{ A'{{1 + eK^Odxu} , dxv). 



Estimate ( 5.4 1 is now a straightforward consequence of Cauchy-Schwarz inequality, and Lemma A.l 



As for the commutator estimates, one uses 

( [A^T]w , v)=eKi{ [A^C]M , V )+iyfien2{ [A'X]dxU , dxV ). 
Estimates (15. 5|) follow, using again Cauchy-Schwarz inequality, and Lemma lA. 21 
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The following lemma offers an important invertibility result on T. 



Lemma 5.3. Let ( £ L°°(R) be such that (H2) is satisfied. Then the operator 

T[eC] : H^{R) — > L^i^) 
is one-to-one and onto. Moreover, one has the following estimates: 
(i) (T[eC])^"'" '■ L^ ^ ^u(If^) ^■s continuous. More precisely, one has 

II 'J" l|L2(R)^ifi(R) < Co, 
with Co = C'(/^o2^<5min''^max,e|CLoo)• 

(ii) Additionally, if C ^ H''°^^{M.) with Sq > ^ then one has for any < s < Sq + 1, 

II '^ llH'^(K)-i.ff;;+^(R) - Cso + l- 

(Hi) If C ^ H'^{M.) with s > So + 1, So > |, then one has 

II '^ llif=^(R)^H;+i(R) - '^s 

where c^ = C(/iq2 , (5^^1,1, (^max, e|C|ff^)' ^^'"■^ uniformly with respect to {iJ,,e,S,j) G VcH- 

Proof. To show the invertibility of T we use the Lax-Milgram theorem. From the previous Lemma, 
we know that the bilinear form: 

a{u, f ) = ( Tu , w ) = ( (1 + eKi()u , u ) + /ii/( (1 + eK2C)dxU , d^v ) 



is continuous and uniformly coercive on (_ff^(M)), provided ( |H2| holds. For any /i > 0, the dual of 
ifHM) is -ff^^(M), of whom L^(IR) is a subspace, and one has {f,g) < |/| rri |<?| r25 independent of 

/d > 0. Therefore, using Lax-Milgram lemma, for all / e L'^{R), there exists a unique u € -ff^(K) 
such that, for all v e -ff/,(K) 

equivalently, there is a unique variational solution to the equation 

(5.6) 1u = f. 
We then get from the definition of T that 

(5.7) vfi (1 + eK2C) dlu = (1 + eKiC)w - fj.ei'K2{dxC){dxu) - f. 



Now, using condition (|H2|, since u e H^{R), ( <^ L°°{R) and / G i^(M), we deduce that d^u G 



L2(M), and thus u G i7^). We proved that T[eC] : H^(R) — > i (K) is one-to-one and onto. 
Let us now turn to the proof of estimates in (i) — (Hi). 



We start from the equality a{u,u) = {f,u). Using elliptic inequality (5.3) and Cauchy-Schwarz 
inequality, one has 

— |"|ffi < a{u,u) = { f , u) < \f\^2\u\^2 < |/L2|w|jyi- 

Co ^ M 

Dividing by Cq u ^ yields the estimate in (i). 

Let us now assume that / e H^{M.), for s > 0. We apply A" to equation (5.6| and we write it 
under the form: 

2•(A''w) = A^/-[A^T]u. 
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Proceeding as above, we use the L^ inner product with A^u, and deduce 
1 



Co 



lA^ul^i < aik'n^K'n) = {%K'u,K'u) = (A^ - [A^2:]u , A^'w) 



(5.8) 



< \A'f\^AA'u 



L2r- ^\L-2 



|([A^X]^i , K'u)\ 



The result is now a consequence of (5.5) 
— //0<s<so + li one has 

1 

Co 
thus 



"iL+i < l/ln^hlff^ + «^o|cla.o+ihlH>lfff+i 



Iffj+i < Co ( |/|^, + eC7o|C|j^so+i|u|^. 



The estimate of (ii) for < s < 1 follows, using estimate (i) and I A'* ^u\f,i ^ l^lffi- ^^^ result 
for greater values of s, 1 < s < sg + 1; follows by continuous induction. 



w|L+i < |/|ff.hlff= + e^o (Id^.o+ihliff + |c|ff=hlff=o+0hlif2+i' 



— //s > So + 1, then plugging ( |5.5| into ( |5.8[ ) yields 
1 

Co 
thus 

I"Ihj+i < Co ( |/|^. + e^o (|C|ff.o+ih|ff, +|C|^.|u|^=o^ 

As above, the result follows by continuous induction on s. 



U 



Finally, let us introduce the following technical estimate, which is used several times in the 
subsequent sections. 



Corollary 5.4. Let C, G W^iM) with s > so + 1; so > 5- Assume moreover that u G H"^ 
that V G H^{M.). Then one has 

(5.9) 



and 



< e C(/lo2\'5^L,^max, |C|fxjkK..-il«l 



[A\1[eC]]l-'H]u , V 

Proof. The first identity can be obtained through simple calculation: using that T[eC] is symmetric, 

{[A^,%-'[eC]]u,T[eC]v) = {lH][^^'^-'H]]u,v) 

= (T[eC]A'"I-i[eC]u - A^w , v) 

= (-[A^T[6C]]T-MeC]«, v)- 



The estimate is now a direct application of ( |5.5[ ) and Lemma 5.3 From point ('zij and (Hi) of 
Lemma |5.3[ one has 



with C = C(/iq2 7 (^niin' '^maxj e|C| tjs-ii £|C| cfsQ+i)- Now apply commutator estimate (5.5|, and one 
obtains straightforwardly our desired estimate. D 



6 Linear analysis 



Let us recall the system ( |4.13[ ) introduced in section 4.2 

= 0, 



(6.1) 



dtC + dx I - — ^. V 



^hi +7/12 
T[eC] {dtv + e<;vdxv) + (7 + S)qi{eC)dxC 



-rh^ 1^12 _ ^1^12 



l+7'»2) 



= -t^elj^d4{dxv)'), 
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with hi = l-e( ,h2 = l/5 + eC, qi{X) = 1 + k^X {i = 1,2) , k,,<; defined in(|4^,(|432]), and 

1[eC]V = qi{eC)V ~ fii^d^ {q2{eC)d.V) , 
In order to ease the reading, we define the function 



One can easily check that 
Additionally, let us denote 



f-x 



hih2 
hi + 7/12 



{l-X)i6-'+X) 
l-X + j{d-^+X)' 



and 



no 



(/ii+ 7/12)2 ■ 



2 1-7 



1, h\-^hi 



so that one can rewrite (here and in the following, we omit the bar on v for the sake of readabihty) 
dtC + fiOd^v + ed^CfiOv = 0, 



(6.2) 



T 



[dtv + ^<;d^{v^)) + {j + S)qi{eC)d^C + eqiieOd^iqsieCy) + fieKd4{d^v)' 



0. 



The equations can be written after applying T ^ to the second equation in (6.2) as 
(6.3) dtU + Ao[U]d^U + Ai[U]d^U = 0, 

with 

where Qo{e(^),Qi{e(^,v) are defined as 

(6.5) Qo«) = (7 + '^)9i«), Oi(eC,«) = '7i«)g3«y 
and the operator n[eC,w] defined by 

(6.6) QH,v]f = 2qi{eC)qs[<)v! + tind^ifd^v). 

Following the classical theory of hyperbolic systems, the well-posedness of the initial value 
problem of the above system will rely on a precise study of the properties, and in particular energy 
estimates, for the linearized system around some reference state U_ = (C, w)^: 



(6.7) 



dtU + Afi[U]d:^U + Ai[U]d:,U = 0; 



U^ 



|t = 



C/o- 



In the following subsection, we construct the natural energy space for our problem. Energy esti- 



mates are then proved in section 6.2 Finally, we state the well-posedness of the linear system (6.7) 
in section l631 



6.1 Energy space 

Let us first remark that by construction, one has a pseudo-symmetrizer of the system, given by 



(6.8) 



'Qo(eC) 







S\U]=\ /(<) 

V ^[<]. 
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Notice that S[U] is symmetric, as T[eQ is symmetric, and one can easily check that S'[[/]Ao[C/] 
is symmetric as well. On the contrary, the operator 



s[u]Mu] 











e^Qi{eC,v) e,1[eC]{v-) 



represents the defect of symmetry. However, T[eC](w-) has the desired following property: 
{l[eC](vd^V),V) = ( qiieQvd^V + f,d,{q2{eC)d^{vd^V)) , V ) 



(6.9) 



2 

UdAqiKhW, V) - fi{q2{eC){d.v)d..v, d.,v) + fil{d,{q2{eCh)d,v, d^v). 



Therefore, the inner product (T[eC,]{vdxV),V) is controlled by fl^ Lri , which is bounded in our 

energy space, as defined below. In the same way, one can control the contribution of e^(3i(eC,w), 
using the smallncss of e through the assumption of the Camassa-Holm regime e = 0{,JJ1). 

Remark 6.1. In the analysis below, the only place where the smallness assumption of the Camassa- 
Holm regime, e = 0{y/jL), is used (apart from the simplifications it offers when constructing sys- 
tem [Q.l]), stands in the estimation of the contribution of e^Qi{e(^,v). As a matter of fact, this 



assumption is actually not required: one could replace Qo(eC) by QoieC,) +^'^Qi{^CtV) *^ t^^ pseudo- 
symmetrizer S[U], thus canceling out the e^Qi{e(^, v) term. Note however that the energy would then 
be slightly different than the one defined below, and that in order for Lemma\6.3\ to hold, one should 



add an additional smallness assumption on ev, in order to ensure Qol^C) + ^ Qii^CiV) > Cq > 0. 

Let us now define our energy space. 

Definition 6.2. For all s > and T > 0, we denote by X" the vector space H^iR) x H''+'^{R) 
endowed with the norm 

yu = ic,v)e x^ |t/||. ^ Id?,. + \v\l. + t^\dxv\jj., 

while X^ stands for the space of U ^ {(, v) such that U E C°{[0, 7]; X") and dtC G L°°{[0, 7] x R), 
endowed with the canonical norm 



\\U\\x^^ sup |C/(i,-)|A- 

te[0,T/e] 



ess sup \dt({t,x) 

te[0,T/e],a:GR 



A natural energy for the initial value problem (|6.7[) is given by 



(6.10) 



E'^iUf = iA'U,S[U]A'U) = (A^C, 



MO 



AX) + A'v 



n<] 



A'v 



/(<) "^' ' \" ^'1 + 6 
The link between E''{U) and the AT^-norm is investigated in the following Lemma. 



Lemma 6.3. Let s > and ( £ L°°{R), satisfying (HI) and (H2|. Then E^{U) is equivalent 
to the I • \x''-norm uniformly with respect to p — (/i, e, (5, 7) G VcH- More precisely, there exists 

Co = C(/loi\ft.o2\'5niax,<5min) > SUch that 



1 

Co 



E'{U) < IC/L, < coE'iU) 



Proof. This is a straightforward application of Lemma 

Qo{<) 



5.2 



(6.11) 



xeR /(eC) 



>(infgo(eC))(sup/(eC))"' 

xes. - xeR 



sup 

xeR 



and that for Qo(eC)i/(eC) > 0, 

MO 



/«) 



<(supgo(eC))(inf/(eC))"\ 

xeR - xeR - 



where we recall that if (HI I is satisfied, then hi ^ 1 — ed^ and h2 



e( satisfy 



inf hi > hoi, 

x€:R 



sup \hi\ < 1 + l/S, inf h2 > ft-oi, sup Ih^l < 1 + 1/(5. 

£ceR ^eR xeR 
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We conclude this section by proving some general estimates concerning our new operators, which 
will be useful in the following subsection. 



Lemma 6.4. Let {ii,e,6,j) — p ^ Vch, satisfying (4.7). Let U = {(u,u)^ such that (u G L' 



< C\V\^„\W\^„, 



satisfies (|H1[),(|H2[). Then for any V,W e X°, one has 

(6.12) ( S[U]V , W 

with C = C{hQ^,hQ2,Ke.^,5:^l^,e\Cu\^^) . 

Moreover, if Cu G H",}^ € X*~^ with s > sq + 1, sq > 1/2, then one has 

(6.13) \[[A^,S[U]]V ,W 

(6.14) ( [A',S~^[U]]V , S[U]W 

with C = C(/loi\ /lo2\ <5max, ^^^j^n' *=|C«|//a) • 

Proof Let U = (Cw,"")^ e X" , V = (C^,^)^ G X^ VK = [Cw^w^ e X". Then by definition of 5'[-] 



< C V 



x^ 



-iFU-o 



< c \v\^^,_,^^^_,\w\^„ 



in (6.8), one has 



( S[U]V ,w) = [ ^^C. , C. ) + ( n<u]v ,w). 



f{<u 

The first term is straightforwardly estimated by Cauchy-Schwarz inequality: 






< 



/««) 



. I Ck I ^2 I C 



U) |2^2 ; 



and t,,^\ is uniformly bounded since Qu satisfies (IHll) 



/«,.) 



The second term is estimated by Lemma 5.2 (5.2). 



|(TK«]t^ , u')! < co\v\mJw\Hi^ <co\V\xo\W\xo. 



Estimate (6.12) is proved 



Now, let us decompose 

( [K',S[U]]V , W 



[k' 



/««) ■ 



C C») + ( [A^T[eCu]]«, w). 



By Cauchy-Schwarz inequality and Lemma |A.2[ one has 






/(eC«) J l-H"'"'l^''l-H"""'l^"'li^ 



where we used Lemma A.l and continuous Sobolev embedding for the last inequality. 



The second term is estimated using Lemma 5.2 (5.5) 



|( [A^T[eC„]]^' , w)\ < C{e\Cu\H^)\V\xs-AW\xo. 



Estimate (6.13) is proved 



Finally, one has 
( [A^5-Mf/]]V^ , S[U]W ) < ( [A^ ^gy]C. , T^C. ) + ( [A^T-M<«]]- , TK«]^^ ). 
The first term is estimated exactly as above, noticing that both /(eCu) a-nd Qo(eC«) £ir6 bounded 



from above and below since (HI) and (H2) are satisfied 



The second term is estimated using Corollary 5.4 (5.9). Estimate (6.14) follows, and the Lemma 
is proved. D 
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6.2 Energy estimates 

Our aim is to establish a priori energy estimates concerning our linear system. In order to be 
able to use the linear analysis to both the well-posedness and stability of the nonlinear system, we 
consider the following modified system 

.g ^gx I dtU + Ao[U]d.U + Ai[U]d^U = F; 

where we added a right-hand-side F, whose properties will be precised in the following Lemmas. 

We begin by asserting a basic X° energy estimate, that we extend to X* space (s > 3/2) later 
on. 

Lemma 6.5 (X° energy estimate) . LetT > andU & L°°{[Q,T/e]\X°) andU,da,Ue L °°{[0, T/e]x 
M) such that dtC € L°°([0,r/e] x M) and ( satisfies ^^,^^, and U,U satisfy system ( |6.15[ ), with 
a right hand side, F, such that 

{F,S[U]U) < CFe\U\l, + fit) |C/|^„, 

with Cp a constant and f a positive integrable function on [0, T/e] . 

Then there exists A = Ci\\dtC\\L^(^[o^T/e]xmr\\0\L^ilo^T/e]xRrPML^i[o,T/e]xMy^F) ««c/i 
that 

(6.16) Vi e [0, -], E°{U){t) < e'^'E°{Uo) + f e'^^'~*'^ f{t')dt'. 

^ Jo 



The constant A is independent of p = {fi^e,S,^) G Vch, o,s defined in (4.7|, but depends on 

Proof. Let us take the inner product of ( |6.15[ ) by S[lf\U: 

{dtU,S[U]U) + {Ao[U]d.U, S[U]U) + {Ai[U]d^U,S[U]U) = {F,S[U]U) , 
From the symmetry property of S[U_], and using the definition of E''{U), one deduces 

(6.17) \j^E\U) = \{U, [duS[US\U) - {S[U]Ao[U]d.U,U) - (5[t/]^i[L/]9,t/,l/) -t- {F,S[U]U). 
Let us first estimate {^S\U]Af)[U]dxU , C/). Let us recall that 

so that 

{S[U]A^[U]d.U,U) =-i (c,e54^^/'(eC)^)C)-(C,a.(Qo(<)>) + e(il[eC,v]d.v,v 

= Ai + A2 + A'i 

The estimates concerning Ai and A2 are straightforward. Using Cauchy-Schwarz inequality, 
there exists C = C{\\U\\ ,^ + \\d,j.U\\ , ^) such that 

IA1H-IA2I < eCdd'. + IH'j < eC|f/|^„. 
As for A3, one has (recalling the definition of £3 in ( |6.6| ): 

{Q[eC,v\d^v,v) = -{dx{qi{eCjq:i{eCjv)v,v) - nK{{d^v){dxv) , d^v) 
Those two terms are controlled, again thanks to Cauchy-Schwarz inequality, so that 

\{£lH,v]d,v,v)\ <C{\\u\\^^ + \\d,U\\^^)\v\\, < C(||C/||^^ + ||9,C/||^^)|f/|^o. 
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Altogether, we proved 
(6.18) 



{S[U]MU]dxU,U)\ < eCi\\U\\^^ + \\d,U\\^^)\U\xo- 



Let us now estimate {S[U]Ai[U]dxU, U). One has 

{S[U]Ai[U]d,U,U) = {e^QM,v)dx(:,v)+e<;{'I[eC]{vdxv),v) 



A4 + A 



4 -1-^5- 



In order to control the term Aa, we write 



{e^QiieCv)d,C,v) = e\qiieQ4{eC)v'd,C,v 



'idAqM)q'M)v'v)X) 



= -e'{d^{qi{eC)q3K)v')v,(:) - e'iqiieQqsKh'd^vX). 



Since p e Vch, as defined in (4.71, one has e < M^fji, and therefore 



\M<^c{\\u\\^^ + \\d,u\\^j\u\^,. 

(where we used, once again, Cauchy-Schwarz inequality.) 

In order to control A^ one makes use of the identity given in ( |6.9[ ), applied to V^ = w, and deduce 
easily 



Altogether, one has 
(6.19) 



lAI < eC-dlC/IL + \\dMr^)E\uf. 



\{S[U]MU]d.U,U)\ < eC{\\u\\^^ + \\d^U\\^^)\u\\,. 



The last term to estimate is {U, [9t, 5[C/]]f/). 
One has 



{U,[dt,S[U]]U) = {v, 9t,S«) + (C, 



dt 



Qo(eC) 



/«) 



C) 



QoiO 



(v,{dtqM))v) + ^l(y,^.,{{^tq2{eO){^,v))) + (c,a,(^^^)C 



e^i 



[v, {dtQvj - fieK2 (dxv, {dtQdxv) 

Qo«)/K)-Oo(eC)/'«) 



+ e 



{<■ 



f{<y 



(dtOC 



From Cauchy-Schwarz inequality and since C satisfies (|H1|), one deduces 
(6.20) 



^{U,[dt,S[U]]U) I < eC{\\dtC\ 



L°=([0,T/e]xR)' llillL°=([0,T/£]xK)''l Ix" 



One can now conclude with the proof of the energy estimate. Plugging (6.18), (6.19) and (6.20) 



into (6.17), invoking Lemma 6.3 and making use of the assumption of the Lemma on F, yields 



1 d 
2~dt 



E'^iUf < eCoE'^{Uf + f{t)E%U), 



where Cq _ C^(|ptC||ioo([o^T/c]xB)' II — llL~([0,T/e]xE)' Ir^— llL~([0,T/e]xI 



,Cf). Consequently 



dt 



E%U) < CoeE'{U) + f{t). 
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Making use of the usual trick, we compute for any A G M, 



at 
Thanks to the above inequality, one can choose A — Co, so that for all t G [0, — ], one deduces 

at 
Integrating this differential inequality yields 



(6.21) 



T 



Vt G [0, -], E%U){t) < e'^'E°{Uo) + / e'^^'-''^ f{t')dt' 



This proves the energy estimate (6.16). 



n 



Let us now turn to the a priori energy estimate in "large" X" norm. 

Lemma 6.6 [X^ energy estimate). Let so > 1/2, s > sq + 1- Let U — {CtV)^ and U_ — {C,v)^ be 
such that U,Ue L°°([0,T/e];X%dtC G L°°{[0,T/e] x E) and C satisfies (|H1]) , (|H2]) uniformly on 
[0,r/e], and satisfying system (6.15) with a right hand side, F, such that 



where Cp is a constant and f is an integrahle function on [0, T/e] . 

Then there exists A — C^MJVi-^^ + Cp) such that the following energy estimate holds: 



(6.22) 



E'{U){t) < e'^'E'{Uo) + / e'^'^'-''\f{t')dt' 

Jo 



The constant A is independent of p = {fi^e,S,"f) G VcH; o,s defined in (4.7), but depends on 

"min'^max, "oi ' ^^02 • 

Remark 6.7. In this Lemma, and in the proof below, the norm \\U\\x-s is to be understood as 

essential sup: 

||[/||js:^ = esssup |C/(i, •)|x= + esssup \dtC{t,x)\. 
te[o,T/e] tela,T/e],xeM. 



Proof. Let us multiply the system (6.15) on the right by A^S[U]A^U, and integrate by parts. One 
obtains 



(6.23) {A''dtU,S[U]A''U) + {A''Ao[U]d^U,S[U]A'U) + {A'Ai[U]d^U,S[U]A''U) 

= e{A'F,S[U]A''U) , 

from which we deduce, using the symmetry property of S[U_], as well as the definition of E''{U): 

(6.24) 
\j^E'{U) - ^(A^C/, [duS[m]A'U) - {S[U]Ao[U]d.A'U,A'U) - {S[U]Ai[U]d.A'U,A'U) 

- {[A\AQ[U] + Ai[U]\d^U,S[U]A'U)+t{A'F,S[U]A'U). 

We now estimate each of the different components of the r.h.s of the above identity. 

• Estimate of {S [U]AQ[U]d.^A''U, A'U) and (5[C/]^i[f/]9^A"[/, A^C/) . One can use the L^ estimate 
derived in (6.18), apphed to A'^U. One deduces 



(6.25) 



\{S[U]A^[U]d.A'U,A'U)\ < eC{\\U\\^^ + \\d,U\\^^)\U\ 
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Now, thanks to Sobolev embedding, one has for s > sq + 1, sq > 1/2 

ci\\u\\^^ + \\d,u\\^j < C(||C/||^J, 

so that 



(6.26) 



\{S[U]Ao[U]d.f^'U,A'U)\<eC{\\U\\^^)\U 



jCrr I \ J\. ^ 



Similarly, using (6.19[), applied to A^U and continuous Sobolev embedding, one has 



(6.27) 



\{S[U]MU]dxA'U,A'U)\<eC{\\U\\^^)\U 






• Estimate of ( [A'' , A [U]] d^ U,S[IJ]]A''U), where A[U] = Aq[U] + Ai[lJ]. Using the definition of 
A[-] and S[-] in |6^,(|6]8|, one has 

{[A^,A[U]\d^U,S[UWU) = ([A^e/'(6C)«]9,C+ [A^/(eC)]9.^' , ^J^A'C 

+ [W,T\Q{<,v)d.,C) , %A'v) + e([A^T-lO[eC,^l]+?^;]a,^;,SA^^;) 

= -Bi + -B2 + -B3. 

Here and in the following, we denote T = T[eC] and Q[eC,,v) ~ Qo(eC) + ^^QiiK^v)- Let us treat 
each of these terms separately. 

- Control of Bi = ([A^e/'(eC)^^]a,C + [A^ /(eC)]a,i; , ^^A%). 

From Cauchy-Schwarz inequality, one has 



\Bi\ < 



[A^e/'(eC)z;]a,C+[A^/(eC)]5.« 



L2 



Q{<,v) 



/(C) 



A'C 



L2 



[A^6/'(eC)^]a.C+ [A^/(6C)]a,z; ^^ < (|9.(e/'(eC))lH=-i + m.fi<))\H^-^) \d.U\H^ 



Since s > Sq + Ij ^^ can use the commutator estimate Lemma |A.2| to get 

<eC{\\U\\^^)\U\^,. 
where we used, for the last inequality. 



It follows, using that r/^'-. € L°° since ( satisfies (HI) 



/(C) 



\Bi\ < eC(||[/|Lj|C/ 



J\.rj^ ' \ I ^ ^ 



- Control of B2 = ([A' ,%-\Q{eC,v)-)]d^(: , "lA'v). 
By symmetry of T, one has 

B2 = (5[A^T-l(Q(eC,«)•)]5xC , A^^^). 

Now, one can check that, by definition of the commutator, 

T[A^T-l(g(eC,«)•)]9.C = %A'T'Q{eC,v)d,C - Q(eC,t^)A^9,C 

= 'IA'T'Q{eC,v)d^C - A'm-\Q{eC,v)d^O + A'{Q{eC,v)d,C) - Q{eC,v)A'd^C 
= -[A^T]T-l(0(eC,«)9.C) + [A\Q{eC,v)]d^C 
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We can now use Corollary |5.4[ and deduce 

The last inequality is obtained using Lemma |A.1| 



Hf, 



From Lemma 



A. 2 and the explicit definition of Q = Qo + e Qi in (6.5), one has 



so that we finally get 

IS2I < eC(||t/||^j|C/|^,. 



- Control of B3 = e('[A^X"^0[eC, w] + ^v]d^v,'IA'vj 
Let us first use the definition of £l[e(,v] (6.5 1 to expand: 



+ e?([A^^;]a,^;,XA^^;) 
= B31 + B32 + -B33. 
In order to estimate -831 and i?32, one proceeds as for the control of i?2- One can check 

%.W .T\qi{<)qM)v-)]d^.v = -[A\%%-\qM)qM)^d.v) + [A^q^{eOqM)v]d,v. 
As above, using Cauchy-Schwarz inequality, Corollary |5.4| and Lemma |A.2[ one obtains 

IB31I < eC(||[/||^j|C/|i,. 
In the same way, 

%[A%T\d^{{d^v) ■)]d^v - -[A',%]%-\d,{{d,v){d.,v))] + d,{[A',d^v]d^v). 

Again, using Cauchy-Schwarz inequality. Corollary |5.4| and Lemma |A.2[ one has 

W,'m-\d^{{d^v){d,v)%A'v) < Csfi\C\Hs\dA{dxv){dM)\Hs-i\y\H'+^ 



— Cs\C\ttbWttS + i\v\ttS + 1, 



and 



A^ 



{d4[A',d,v]d,v),A'v) = fi{[A',d,v]d,v,A'd^v 



— Cs\v\ ttB + 1 \v\ „s + 1 



Thus we proved 



IB32I < eC(||t/|Lj|C/ 



5.2 



one has 



Finally, we turn to i?33 = ec^l [A'^,y]dxV,'I_A^v]. From Lemma 

1^331 < \[A\v]dxv\^,\A%\^, < \[A\v]dxv\^,\A'v\^, + y/Jl\dx{[A\v]dxv)\^,\A%\^,. 
Note the identity 



so that Lemma 



A.2 



94[A^^;]9,l;) = {[A' , d^vid^v) + {[A' , v]dlv) , 
yields y^lOxdA" ,v]dxv)\^2 < |:y|^=+i|w|^=.+i- 



Altogether, we proved 



(6.28) 



\{[A',A[U]]dxU,S[U]]A'U)\<eC{\\U\\^/)\U 
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• Estimate of ^{A^U, [dt,S[U]\A'U). 
One has 

(A^C/, [dt,S[U]]A'U) = {A'v, IdtMA^v) + (A^C, 



Qo(<) 



A^C) 



= A 



e^i 



'«, (9t(?i(eC))A^t;) + ^i(^A'v,d4{dtq2{eC)){d,A'v)) 
{^A'v, {dtCWv) - fien^i^A'd^v, (dtQA^d^v)) 
+ e(A^C, = -^^^^p ^{dtCWC 



(A'C,ft(^)A-< 



From Cauchy-Schwarz inequality and since (^ satisfies (HI I, one deduces 
, 1 



-{A'U,[dt,S[U]\A'U) I < eCiptC 



C 



and continuous Sobolev embedding yields 

, 1 



(6.29) 



{A'U,[dt,S[U]]A'U) \ < eC(||C/||^J|C/ 






U\ 



!lllL°°([0,T/e]xR)' \\^\\L°°([n,T/e]xR)'\ \X 



Finally, let us recall the assumption of the Lemma: 



(6.30) 



{A'F,S[U]A'U) <eCF\U\l^, + f{t)\U\ 



js:= 



We now plug (6.261, (6.27), (6.281, (6.291, and (6.301 into (6.24 1. Using Lemma 6.3, it follows 



2di 



E'iUy < CQeE'{Uy + E'{U)f{t), 



with Co — C{\\U_\\ + Cp), and consequently 



-£;■'(£/) < eCoE^iU)+m 



Now, for any A G M, one has 



e'^*dt{e-'^'E'{U)) = -eXE'{U) + -j-E'iU) 



Thus with A ~ Co, one has for all t G [0, — ], 

d 



dt 



{e-'^'E\U)) < J{t)e 



dt 



e\t 



Integrating this differential inequality yields 



T, 



Vi € [0, -], E^U){t) < e'^'E'iUo) + Co I e'^<-'-''^ f{t')dt' . 



This concludes the proof of Lemma 6.6 



n 
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6.3 Well-posedness of the linear system 

Let us now state the main resuh of this section. 



Proposition 6.8. Let sq > 1/2, s > sq + 1? o-nd p = {fi,e,S,j) G Vch, defined in (4.7|. Let 
also Vy — {C^ tV^)^ G X^, as defined in Definition 6.2 he such that (|H1[),(H2 1 are satisfied for 



t S [OjT/e], uniformly with respect to p £ VcH- For any Uq G X"^, there exists a unique solution 
to ^, C/P e C°{[0,T/e];X'')r]C\[0,T/e];X^-^) C X^, and Cq^Xt = C(||C/||^J, independent 
of p (z Vch, such that one has the energy estimates 



V 0<i< 



T 



£^'(C/P)(t) < e^^^*£:''(C/o) and E'-'iOtW^) < Coe'^^*E'{UQ) 



Proof. Existence and uniqueness of a solution to the initial value problem (6.7) follows, by standard 
techniques, from the a priori estimate (|6.22 1 in Lemma 6.6 



(6.31) 



E'{U){t) <e''^'E'iUo), 



(since _F = 0, and omitting the index p for the sake of simplicity.) We briefly recall the argument 
below. 



First, let us notice that using the system of equation (6.7), one can deduce an energy estimate 
on the time-derivative of the solution. Indeed, one has 

\dtU\^,_, = \ - Ao[U]d.U - Ai[U]d^U\^^_, 
< \ef{eQvd:,v + f{eOv^,v\f,.-^ 



(6.32) 



< Ci UL.MUL., < C„e'^^'E'{Uo), 



where we use Lemmata lA.il and 15.31 

The completion of the proof follows by classical methods, that we briefly recall. For more 



details, we refer to [32[|37| , for example. In order to construct a solution to (6.7), we use a sequence 
of Friedrichs moUifiers, defined by J^ = (1 — vd^)~^^^ [y > 0), in order to reduce our system 
to ordinary differential equation systems on X^ , which are solved uniquely by Cauchy-Lipschitz 



theorem. Estimates ( |Ol] ) , ( [6321 ) hold for each U^ G C°([0,r/e];X"), uniformly in t/ > 0. One 
deduces that a subsequence converges towards U G L^([0, T/eJjX**), a (weak) solution of the Cauchy 
problem (6.7). By regularizing the initial data as well, one can show that the system induces a 



smoothing effect in time, and that the solution U G C°([0, T/e]; X^) n C^([0 T/e]; X**-!) is actually 
a strong solution. The uniqueness is a straightforward consequence of (6.31) (with Uq = 0) applied 
to the difference of two solutions. D 



7 Proof of existence, stability and convergence 

In this section we prove the main results of this paper. We start by proving an a priori estimate 
on the difference of two possible solutions. Existence and uniqueness of the solution of the Cauchy 
problem for our new Green-Naghdi system, in the Camassa-Holm regime e — 0{y/Ji) and over large 
times, is then deduced from the linear analysis of the previous section and this a priori estimate. 
The estimate also provides 

• the stability of the solution with respect to the initial data, thus the Cauchy problem for our 
system is well-posed in the sense of Hadamard, in Sobolev spaces; see subsection |7.2[ 



the stability of the solution with respect to a perturbation of the equation, which allows, 
together with the well-posedness, to fully justify our system (and any other well-posed, con- 
sistent model); see section [s] 



April 17, 2013 



Vincent Duchcnc, Samcr Israwi, Raafat Talhouk 



7.1 One more a priori estimate 

In this subsection, we control the difference of two solutions of the nonlinear system, with different 
initial data and right-hand sides. More precisely, we prove the following result. 

Proposition 7.1. Let sq > 1/2, s > Sq + 1, and assume that there exists Ui for i € {1,2}, such that 
t/, = (Cw.)^ e X^, C/2 e L°^{[0,T/e);X''+^), Ci satisfies ^,^ on [0,T/e], with /ioi,/io2 > 0, 
and Ui satisfy 

dtU, + Ao[U,]d,Ui + A,[U,]d,U, = F, , 
dtU2 + Ao[U2]d^U2 + Ai[U2]d^U2 = F2 , 

with Fi € L°°{[0,T/e];X'^). Then there exists Co = C{h^^ , /ig2 , i^max, <5,7iin, e|C^i| j^-s, e|C^2| -v-s) and 
At = Co X C(|C/2|l~([o,t/£);X=+i)) such that 

Vte[0,-], E%Ui-U2){t)<e'^^'E%Ui\t=o~U2\t=o) + Co f e'^^^'-*'^E%F, ^ F2){t')dt' . 



Proof. When multiplying the equations satisfied by Ui on the left by S[Ui], one obtains 

S[Ui]dtUi + Eo[[/i]5,[/i + Si[C/i]9,C/i = S[Ui]Fi 
S[U2]dtU2 + J:o[U2]d:,U2 + Si[t/2]a,C/i - S[U2]F2; 

with So[C^] — S[U]Ao[U] and Si[[/] — S[U]Ai[U]. Subtracting the two equations above, and 
defining V — Ui ~ U2 = (C, v)^ one obtains 

s[Ui]dtV + So[c/i]a,y + ^i[Ui]d^v = 

S[Ui]{Fi - F2) - (I]o[(7i] + I]i[C/i] - So[C/2] - ^i[U2])d^U2 - {S[Ui] - S\U2]){dtU2 - F2). 



We then apply 5 ^ [Ui] and deduce the following system satisfied by V: 

(7.1) 

where, 



dtV + Ao[Ui]d^V + Ai[Ui]d^V = F 
V{0) = (Ui - c/2) |t=o , 



(7.2) F = Fi-F2 - 5-i[[/i](Eo[C/i]+Si[[/i]-Eo[C/2]-Ei[C/2])9,f/2 

- S-^[Ui]{S[Ui]~S[U2]){dtU2^F2). 



We wish to use the energy estimate of Lemma 6.6 to the linear system ( |7.1[ ). Thus one needs to 
control accordingly the right hand side F. 

In order to do so, we take advantage of the following Lemma. 



Lemma 7.2. Let s > sq > 1/2, and (/i,e,(5,7) = p € 'Pch, satisfying (4.7|. Let V = {(^v,v)^ , 
W — {(yj,w)^ e X'^ and Ui — (Cij"!^!)^; U2 — (C2,W2)^ £ X^ such that there exists h> with 



l~eCi>h>0, 1 - €(2 > h > 0, ^ + eCi > h > 0, ^ + eCi > h > 0. 

d 



Then one has 
(7.3) 



A' {S[U,] - S[U2])V , W 



< eC \Ui~U2\^,\V\^^\W\^, 



(7.4) i^A'{S[Ui]A[Ui]-S[U2]A[U2])V , W j < e C |C/i - C/2|^,|l^|^, |T4^|^o 

with C = C(/i"\(5max,i5min''^|^i|jf='«|^2|_YJ, and denoting A[-] = Ao[-] + Ai[-]. 
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Proof of Lemma\T^ Let V = {Cv,vy, W = (Ct«,w)^ € X° and Ui = {Ci,viV, U2 = ((2,^2)^ G 
X^ . By definition of 5[-] (see (6.8 1), one has 



(^k'{S[Ui]^S[U2])V , 



w 



A' 



'QoieCi) Qo«2) 



)C. , C») + {A'{l[<i]~n<2])v, w). 



/«l) /«2) 

Now, one can check that 

neCl]v-1H2]v = {qMi)-qM2))v-f^>^d4{q2{eCi)~q2{eC2))d,v} 
= e('«i(Ci -C2)i'-M'^52;{k2(Ci -C2)9xw}j, 



SO that, after one integration by part, and using Cauchy-Schwarz inequality and Lemma A.l 
has 



one 



(7.5) 



( A"(l[eCi]-1[eC2])w , w) < e C(Ki,t/K2)|Ci-C2|^.|w|^=+i|w|j^i. 



In the same way, we remark that one can write f°(x) ^^ 9- rational function: 

QoiX) ^ ^^^^^,{l + X)il-X + jiS-'+X)) _ P{X) 



fix) 



il-X){S-^+X) 



Q{x)- 



It follows 



QoKi) _ Q0K2) ^ PKi)QK2)-P(eC2)QKi) 
/(eCi) /«2) Q«i)QK2) 

^ (P(eCi) - P{eC2))Q{<2) - P{<2){Q{<i) - g(6C2)) 

Q«i)Q«2) 



Since P(A) and Q{X) arc polynomials, and using Lemma [A.l it is straightforward to check that 
for s > sq > 1/2, one has 

|(P(6Ci)-P(eC2))Q(eC2)-P(eC2)(QKi)"Q(eC2))|^. < ^^(elCil^., e|C2|^j|Ci - (2!^. 



Now, applying Cauchy-Schwarz inequality and Lemma |A. 1| together with Lemina |4.5[ one deduces 
that as long as 

hi{eCi),h2{<i)M{<2),h2{<2) > h > 0, 

one has (again thanks to continuous Sobolev embedding for s > sq > 1/2) 
'QoieCi) Qo(eC2)^ 



(7.6) 






•sW \j^2 1 



with C = C(/i-i,e|Ci|^.,e|C2|^.,'5max,'5^i\j. Estimates ^ and ^ yield (|73|. 



Let us now turn to (7.4 1. One has 
(7.7) 
[a^{S[UM[Ui]-S[U2]A[U2])V, w) = e(A^(^^/'«i)«i-^^/'«2)«2)C., C.) 

+ (A^(Qo(eCi)-go(eC2))^', C«.) 

+ (A^(go(eCi)-go(eC2))C., w) 

+ €^(^A'{Q,{eCi,vi)-QM2,V2))Cv, w) 

+ e(^A'{£l[eCi,vi]~£l[eC2,V2])v , w) 

+ e^( A^(X[eCi](fi I') -TK2](^;2 v)) , wy 
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The second and third terms in the right-hand-side of (7.71 may be estimated exactly as in (7.61, 



and we do not detail the precise calculations. The first term follows in the same way, using the 
decomposition 



, ?o(eCl) rl, /- ^ Qo(eC2) fl, ^ ^ \ 

e( ■■- - ~ I (eCi)t^i - -fr^r^f (eC2)i'2) 



/(eCi 



so that one has 



/«2 



''^"^'^^^-/'Ki)-%^/'K2))(..i) 



/«l) 



A ( „. . s / (eCi)wi jny^J (eC2)w2JC« , Qu 



/«1 



f{<2) 



/«2) 

+ <i{vi-V2)-rr---f K2), 

/le'^2J 



- C'l^hllifJ^ICl -C2\hs\C.v\hs\C.w\l2 

+ C{e\C2\^s)(\vi -V2\jj,\C,v\jjs\C,w\^-2- 



The fourth term is similar, as f^Qii^CiT'^i) — QiC^CiiC^i) is a bivariate polynomial. Let us detail 
the last two estimates. One has 

(7.8) {£}[e(i,Vi]-£}[e(2,V2])v = 2{qi{e(i)q3{€Ci)vi-qi{eC2)q3{eC2)v2)v + fiKd^{vd^{vi-V2)). 



Again, the contribution of the first term in (7.8) is estimated as above (recalling that this term is 



multiplied by a e- factor), and the contribution of the last term in (7.8) is estimated below 



e/iK 



(^ A''d^{vd^{vi -V2)) , w ^ 



< Cenlvi — V2\ 



Hi 



\H''\ \m 



We conclude by estimating the last term in (7.7). One has 



T[eCi](wi v) -1[eC2](w2 v) = (gi(eCi)wi ~ qi{eC2)v2)v - ^iv^x{{q2{|^C,l)^x{vl v) - q2{K2)dx{v2 v)] 

= (TKi]-'r[eC2])(i'i v) + e{vi~V2){qi{eC:2)v) 
- tivdx{q2{eC,2)dx{{vi - V2) v)]. 

One finally uses Cauchy-Schwarz inequality. Lemma | A. 1 as well as ( |7.5[ ), and obtain 

li K" {%[eC,i]{vi v) ~ 'Z[eC,2]{v2 v)) , wj < e^ C(ki, i/K2)|Ci - (2!^. |wi w|j:^.+i |w|^i 

-f e C(Kie|C2|j:^j|^'l - ^2|^e|w|jJ. |''«|i2 

-I- ve C(K2e|C2|^j|^'i - "2|jj=+i|w|^=+i|w|^i. 



Altogether, we obtain (7.4), and the Lemma is proved. 



D 



Let us continue the proof of Proposition |7.1[ by estimating F defined in (7.2 1, that we recall: 



(7.9) F = F1-F2 - 5-i[[/i](l]o[C/i]+Si[t/i]-So[C/2]-Si[C/2])a,C/2 

- s-^[u{\{s[u{\- s[U2]){dtU2~ F2). 

More precisely, we want to estimate 

( k'F , SpiWV ) = ( A'^Fi - A''F2 , S[Ui]K'V ) 

+ ( A^(So[f/i] + I]i[C/i] - So[[/2] - ^i[U2])dxU2 , K'V ) 

- ( A^(5[C/i] - S[U2WtU2 - F2) , A^^ ) 

+ ( [A^5-l[^^l]](So[t^l]^-Sl[[/l]-I]o[C/2]-Sl[C/2])5,f/2 , S[U{\K'V) 

- ( [k\S-^[Ui]]{S[Ui] ~ S[U2]){dtU2-F2) , S[UiWV ) 
^{I) + {II) + {III) + {IV) + {V). 
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Let US estimate each of these terms. The contribution of (/) is immediately bounded using 
Lemma 16.41 



(7.10) 



(/) I < C\Fi-F2\xs\V\xs, 



with C^C{h-\Sn,,^,S-l,e\uA. J. 



The contributions of (//) and (///) follow from Lemma 7.2 Indeed, recalling that V = Ui 
U2, ( |7.3[ ) yields immediately 

(7.11) \{III)\ < Ce\dtU2-F2\xAV\'x., 



and (7.4) yields 



(7.12) I (//) I < Ce\d,U2\xAV\l,., 

with C = C^ C(/l-l, (5,nax, Cn. A^X^^A^Ax^)' 



Finally, we control {IV) and {V) using Lemma [6^ (6.14): 



|( [A^5-l[C/l]]f/, S[UA^'V)\<C\U\hs~AV\xs, 

with C - C(/l-l,(5„,ax,Cn,e|Cl|^J- 

Thus it remains to estimate |C/|//=-i, where [/ = U^^) = (Eo[C^i] +i;i[C/i] -So[t/2] -Si[C/2])9^[/2 
or [/ = t/(,,) = {S[UA - S[U2]){dtU2 - F2). 

We proceed as in Lemma [7?2| helped by the fact that one is allowed lose one derivative in our 
estimates. Let W = dtU2 — F2 = {(w,w)^ . One has 



U^u) = {Sm - S[U2])W = 



/( QoieCi) Qo{<2) \ 
V fieCi) f{<2) J^"" 



y (^T[eCi]-TK2]j^ 






Recall that 



T[eCi]w - T[eC2]w = fi^niid - C2)w - fj,vdx{K2{Ci - C2)dxw}j, 
so that one has straightforwardly 

|w(ji)|ff=-i < eC{Ki,iyK2)\Ci -C2|h=|w|^=+i. 



As for the first component, we apply (7.6) and deduce 

IQiolff"-! "^ (^''^ C(u),^^^ C{ii)) < eC'ICi ^ C2|_ff=-i|C«;|//=-i|C(M)lff=-i- 
It follows 
(7.13) I (y) I < Ce\dtU2~F2\xAV\h, 

with C = C(/l~\(5max,5m,?n><^|^l|x='^l^2|^J- 

Now, recall t/(i) = (I]o[C^i] + Si[C/i] — Eo[C^2] — ^i[U2])dxU2- Proceeding as above, one obtains 

I C/w \h.-i < C\dxU2\xAV\xs, 
and thus 



(7.14) 



{IV) I < Ce\dxU2\xAV\ 



X" 



Altogether, we proved (using Lemma 6.3) that F, as defined in ( |7.2[ ), satisfies 
(7.15) \(a'F,S[Ui]A'v)\ < Co{\dxU2\xs + \dtU2 ~ F2\x^)eE%V)^ + CoE%V)E'{Fi ^ F2). 



with Co = C{h ^, (5max, ^min' ^| ^1 1x^2. ^| ^2 1 jf s )• Notice also that by the system satisfied by U2, one 
has (see detailed calculations in (6.32)) 

\dtU2-F2\x^ = ~\{Ao[U2] + Ai[U2])dxU2\x' 
<C(|f/2U=+0 

n 



We can now conclude by Lemma 6.6 and the proof is complete 
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7.2 Well-posedness result 

In this section, we prove the weh-posedness of the Cauchy problem for our new Green-Naghdi 



model (4.13) in the sense of Hadamard. Existence and uniqueness of the solution is given by 



Theorem |7.3[ while the stability with respect to the initial data is provided by Theorem |7.4[ These 
results correspond to Theorems |3.2| and 3.3 as stated in section [3] 



Theorem 7.3 (Existence and uniqueness). Let sq > 1/2, s > so + 1; '^d let Uq — (Coi'^o)^ G X'^ 
satisfy (|H1|,(|H2|. Then there exists a maximal time T^ax > 0. uniformly hounded from, below 



with respect to p £ VcH, such that the system of equations (4.131 admits a unique solution U 



iCvy e CO([0,r„,ax/e);^^)nCi([0,T„,ax/e);^^"') with the initial data (Co,wo)|t=o = (Cc^'o) 



and preserving the conditions (H1),(H2) (with different lower bounds) for any t € [0,Tmax/e)- 

Moreover, for any <T < T^ax, there exists Co, At = C{hi^^ ^h^^ ,Sm!i^,S^^l^,M,T,\Uo\^^), 
independent of p £ Vch, such that one has the energy estimates 



T 

VO<i< - 

e 



At* 



C/(t,.) _ + 5,C/(i, OL.-i <Coe'^- 



If Tniax < oo, one has 



\U{t,-)b 



oo as t 



T 



or one of the two conditions ( |H1| ,( [H2| ceases to be true as t — > T^nax/^- 

Proof. We construct a sequence of approximate solution ([/" = (C"i'""))n>o through the induction 
relation 



(7.16) 



[/" 



Uq, and Vn e N, 



\t = Q 



A[u-^]d,u 



n+1 _ 



0; 



f/o. 



By Prop osition [64 there is a unique solution [/"+i e C°([0, r„+i/e];X") n Ci([0,r„+i/e];X"-i) 
to ([7l6| if [/" e C°([0,r„/e];X'*)nCi([0,T„/e];X*-i) C X|,, and satisfies pn|),(|H2|). 



Existence and uniform control of the sequence f/". The existence of T' > such that the sequence 



[/" is uniquely defined, controlled in X|,, , and satisfies (H1),(H2), uniformly with respect to n G N, 
is obtained by induction, as follows. 
Proposition |6.8| yields 



(7.17) £:^([/"+i)(t) < e'^-^E'iUo) and \dtU''+^{t,-)\^^_, < CnE\U'^+^) < C„e^^"*£;^(C/o), 



with A„,C„ — C(||C/"||^^ ,hoi,n,ho2,n), provided C/" € X^^ satisfies (H1|,(H2) with ft.oi,n,^02,n 
on[0,T„/e]. 



Since C/" — (C",w") satisfies (7.161, one has 



dtC^' = -e/'(eC")«"5,C 



n+l 



f{eC)d.v 



n+l 



Using continuous Sobolev embedding of i/'* ^ into L°° (s— 1 > 1/2), and since C" satisfies (HI I, (H2) 
with /ioi,n, ^02, n On [0,T„/e], one deduces that 

(7-18) \dtC^'\L^ < C{h^lu, Kin, Cn: '^max) ||C^1Lj^ • 

Let g"+i = a + 6eC"+\ where (a, h) e {(1, -1), (i, 1), (1, ki), (1, ks)}. One has 

Jo 



„n+l 



'■n+l 



SO that ( |7.18[ ) yields 

|5"+i - ff"+M*=o |l~ < ei X 6C(Vi!„,V2!„,Cn,'5-L)||t^"L= 
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Now, one has 5"+^ |t=o = 5" |t=o > niii^(^oi^O' ^02^0) > 0, independent of n. Thus one can easily 
prove (by induction) that it is possible to chose T' > such that g""*"^ > a/2 holds on [0,T'/e], 
and the above energy estimates hold uniformly with respect to n, on [0,T'/e]. 



More precisely, one has that C" satisfies (HI l,(H2) with /ioi/2, /102/2 > and the estimates 

(7.19) E'{U''){t)<e'^'E%Uo) and \dtU^{t,-)\^^_, < Coe'^'E'{Uo), 

on [0,T'/e], where A, Cq — C{\Uo\j^^, h^-^ q, /iq2 0) ^^^ uniform with respect to n. 

Convergence of C/" towards a solution of the nonlinear problem. We now conclude by proving that 
the sequence t/" converges towards a solution of our nonlinear problem. In order to do so, let us 
define V"" = [/"+^ — [/". Vn satisfies the system 



(7.20) 

where, 

(7.21) 

F" = -S-i[C/"](So[C/'^ 



9*1/" + ^o[(7"]a,y + Ai[u^]d^v = f" 

^|*=o =0, 



] + Ei[C/"] - Eo[t/"-i] - El [L/"-i])9^[/" - 5-1 [t/"](S'[C/"] - S'[C/"-i])atf/" 



We wish to use the energy estimate of Lemma |6.5| to the linear system (7.1). Thus one needs to 
control accordingly the right hand side F"^ . 
More precisely, we want to estimate 



( A'^i^" , S[U'']V'' )=-{ (Eo[C/"] + Ei[C/"] - Eo[L/"~^] - Ei[[/"~i])5:,f/" 
- ( (S'[C/"] - 5[[/"-i])at[/" , T/" ). 

Proceding as in Lemma |7.2[ one can easily deduce 



y" ) 



\{k'F'' , S[U'']V'^ )| <eC|t/"-t/" 



1^0 1 



\W^'- 



,e C/" 



Ivi/i. 



J- 



1X0 



a^c/" 



1X0'' 



with C = C(Vl\ ^2': ^max, Cn> e|f^" 

Using the uniform control of [/" , Sj t/" in ( |7.19 1 , one deduces 
with Cq independent of n. Thus Lemma |6.5| yields 



<eCo|F"'i|^o|F"|^„, 



Vi e [0, — ], £;"(y")(i) < eCa I e^^(*-*')S°(F"-i)(t')dt' < eC / £;°(F"-i)(t')di', 

e Jo io 



where C is independent of n and t. Hence 

Vt e [0, — ], £;"(y")(i) < 



e"^"^ 



sup E^V'^W), 

i'G[0,T'/e] 



and the sequence [/" = C/" + ^ F" converges in C"'([0,r7e]; AO). 

Completion of the proof. Since C/" converges in C*'([0,T'/e]; A°) and is uniformly bounded in A*', 
standard interpolation arguments imply that the sequence C/" converges in C°([0, T'/e]; A" ), for any 
s' < s. Similarly, dtU" converges in C°([0,T'/e]; A" ~^). Choosing s' — l > 1/2, one may pass to the 



limit all the terms in system (7.16), and one deduces that the limit [/ is a solution of system (4.131 



Passing to the limit the properties of C/", and in particular energy estimates (7.19), yield U € 

L°°{\0,T/e];X''),dfU G L°°{[0,T/e];X''^) 



and U satisfies the energy estimate of the Theorem 
(using Lemma 6.3), and preserves the conditions (H1),( [H2| ) for any t £ [0, "^ ^ 
P e VcH- 



], independently of 



Finally, as in Proposition 



6.8 



(with f/ = [/), one has f/ e C"([0,rVe]; A'')nCH[0,T7e];A''-i). 
The uniqueness of U follows from the stability result of Proposition |7.1| with Fi = F2 = 0, and one 
can therefore define a maximal time of existence of the solution, that we denote Tmax/f- ^niax is 
bounded from below by T' > 0, and the behavior of the solution as < — > Tmax/e if ^niax < 00 follows 
from standard continuation arguments. D 
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Theorem 7.4 (Stability). Lef s > sq + I wif/i sq > 1/2, andjet Uq,i — (Co,i, ^0,1)^ G -'^'^ a'^c? L/0,2 = 
(Co,2, "^0,2)^ € X^"*"^. Under the assumptions of Theorem 7.3. let Uj be the solution of system (4.13) 
with Uj \t=o = Uqj . Then there exists T, X, Cq — C(/ig2 , h, 
such that 



x= 



|t/o,2|x»+0>0 



yt e [0, 



[Ui - C/2)(t, OL. < Coe'^^'Uifl - C/: 



ix= 



\x= 



Proof. The existence and uniform control of the solution Ui (resp. U2) in L°°([0,T/e]; X*) (resp. 
L°°([0,T/e];X^+i)) is provided by Theorem [73] (choosing T < T^ax)- The proposition is then a 
direct consequence of the a priori estimate of Proposition |7.1[ with Fi ~ F2 — 0, and Lemma [6?3l D 



8 Full justification of asymptotic models 

We conclude our work by explaining how the results of the previous sections allow to fully justify 
our system (and other consistent ones) as asymptotic model for the propagation of internal waves. 
A model is said to be fully justified (using the terminology of [26]) if the Cauchy problem for both 
the full Euler system and the asymptotic model is well-posed for a given class of initial data, and 
over the relevant time scale; and if the solutions with corresponding initial data remain close. As 
described in [28] Section 6.3], the full justification of a system (S) follows from: 

• (Consistency) One proves that families of solutions to the full Euler system, existing and 
controlled over the relevant time scale satisfies the system (S) up to a small residual. 

• (Existence) One proves that families of solutions to the full Euler system as above do exist. 
This difficult step is ensured by Theorem 5 (or Theorem 6 for large times) in [28] , provided 
that a stability criterion is satisfied (see details therein). 

• (Convergence) One proves that the solutions of the full Euler system, and the ones of the 
system (S), with corresponding initial data, remain close over the relevant time scale. 

The last step supposes that the Cauchy problem for the model is well-posed, and is a consequence 
of the stability of its solutions with respect to perturbations of the equation, so that the first two 
steps of the procedure (consistency and existence) yield the conclusion (convergence), and therefore 
the full justification of the model. Let us refer to Theorem 7 in [28] for the application of such 
procedure for the full justification of the so-called shallow-water/shallow-water asymptotic model, 
which corresponds to our system, when withdrawing 0{fi) terms. 



The consistency of our model has been given in Theorem |3.1[ The well-posedness of the Cauchy 
problem is stated in Theorem |3.2[ and the stability results is a consequence of Proposition |7.1[ 



Thus we have all the ingredients for the full justification of our model, stated in Theorem 3.4 and 
that we recall below. 



Theorem 8.1 (Convergence). Letp e Vch (see (i.l)) and C/° = (C°, V'")^ G H'^'^ , N sufficiently 



large, satisfy the hypotheses of Theorem 5 in ^281 as well as (H1),(H2). Then there exists C, T > 0, 
independent o/p, such that 



• There exists a unique solution U = (Ci'^)^ io the full Eu ler system (2.4), defined on [0,T] 
and with initial data (C°,'0*')^ (provided by Theorem 5 in '28^); 



• There exists a unique solution Ua = (Ca,^a)^ to our new model (3.2), defined on [0,T] and 
with initial data (C^i'i'^)^ (provided by Theorem 3.2); 



With V = u[C,?/;], defined as in (3. 3), one has 

|(C,«)-(Ca,t'a)L„([0,T]^j,.)<C(A*2^Bo-l)t. 



Proof. As stated above, the existence of U is provided by Theorem 5 in [28], and the existence of Ua 



is given by our Theorem 3.2 (we choose T as the minimum of the existence time of both solutions; 
it is bounded by below, independently of p G Vch)- If N is large enough, then U = (C, V')^ satisfies 
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the assumptions of our consistency result, Theorem 3.1 and therefore {C,v)^ solves (3.2) up to a 
residual R = (ri,r2)^, with |^| roonn t^i.us-, ^ C{ii^ +Bo^^). The result follows from the stability 



Proposition 7.1 with Fi = (ri,1[eC]"^?-2)^ and F2 = 0. D 



In addition to allowing the complete, full justification of our model, the results of the previous 
sections allow to rigorously justify any lower order, well-posed and consistent model. We quickly 
show how to apply the procedure to such models, with the example of the so-called Constantin- 



Lannes decoupled approximation, introduced by one of the authors in 19 , and that we recall 
below. 

Definition 8.2 (Constantin-Lannes decoupled approximation model). Let C°,ti° be given scalar 
functions, and set parameters {iJ,,€,6,^) G Vch, o.s defined in ( |3.1| ), and (A, 6*) e R"^ . The 
Constantin-Lannes decoupled approximation is then 

UcL = [v+{t,x-t)+v^{t,x + t),{-f + d){v+{t,x-t)-v^{t,x + t))y 

where v± \t=o = |(C° ± ^+3) \t=o o,i^d ^± = (1 i iiXd'^)^^v^ with v\ satisfying 

(8.1) dtvi ± taiv^d^vi ± e'a2{vifd^vi ± e^al'^^ivlfd^Vi 

± ^ivl'^^lvi - ^ivt'^^l^tvi ± iied,{4^\idlvi+4{d.,vif) = 0, 

with parameters defined as follows: 

«! = ;: — T' "2 = -3^ , ex9 ' "3 == -5- 



2 7 + <5 ' " {l + Sy " h + 5) 



i.2) 



"' -6J(<5 + 7) ' "^ " 6 S{S + j) ^' 

,,, ^ {l + j5){d^-j) l~l (1-7) , ,3 (52-7 

1 " 3,5(7 + (5)2 4 ' 6(7 + (5) 27 + 5' 

e ^ (l+7^)(>52-7) .. , 1^. _ (1-7) 

^"2- 3,5(7-^5)2 ^ ^ 4 ^ 12(7 + 5)' 



Remark 8.3. The scalar equation (|8.1f has been introduced as a model for gravity surface wave 
(one layer of homogeneous fluid) in [2^ 1, and its rigorous justification has been developed in UW- 
The justification is to be understood in the following sense: if one chooses carefully the initial data 
(so as to focus on only one direction of propagation), then the Constantin-Lannes equation provides 
a good (unidirectional) approximation of the flow. The reason why a bi-directional, decoupled ap- 
proximation as above has not been developed in the water-wave setting is that lower order scalar 
equations offer the same accuracy. The specificity of internal waves lies in the existence of a critical 
ratio (S^ = 7J for which quadratic nonlinearities vanish, thus calling for higher order decoupled 
models, especially in the Camassa-Holm regime. 



The Cauchy problem for the equation (8.1 ) has been proved to be locally well-posed in [13], and 
a property of persistence of spatial decay at infinity has been proved in flO] . We state these results 
below. 

Proposition 8.4 (Well-posedness and persistence). Let u^ G _ff*+^, with s > sq + 1, sq > 1/2. Let 



the parameters in (8.2) be such that ^,e,z/^' > 0, and define m> such that 



(8.3) ^f'V(//^')-i+A^ + 6+|ai| + |a2| + |a3l + l^f''l + l4'l + l'^2l < m. 

Then there exists T = C{m, \u°\^s+i) and a unique u G C°([0,T/e); iJ^+i) n Ci([0,T/e); iJ^) such 



that u satisfies (8.1 1 and initial data u 



|t=o 



Moreover, u satisfies the following energy estimate for < f < T/e; 

ll^*"llL~([0,T/e);HS) + II"IL-([0,T/c);H;+1) - C'(to, ju^ |^, + i ). 
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Assume additionally that for fixed n, fc G N, the function x^vP G H^^^ , with < j < n and 
s = fc + 1 + 2(n — j). Then there exists T — C(rn,n,k,Y^^=o \x-' u'^\ Tr!!+k+i+2(n-j)) such that for 
< t < T X min(l/e, 1/^), one has 

n 



The decoupled model of Definition 8.2 is consistent with our Green-Naghdi system (3.2), in the 
following way. 

Proposition 8.5 (Consistency). Let(°,v° G il*+^, with s > Sq + I, sq > 1/2. For {^,€,6,j) = p G 
VcH , CLS defined in (|3.1[), we denote Uq^ the unique solution of the CL approximation, as defined 



in Definition 8.2 For some given M*_^_q > 0, sufficiently large, assume that there exists T* > and 
a family ([/^^)pg-p such that 

T* = max{T>0 such that ||t^cL||L~([o,T);ff=+6) + ll^*^CL|Lcx=([o,T);ff=+=) - ^-^^+6 } • 



Then there exists V^ = [/"^[C/p^] such that U — U^j^ + U'^ satisfies the Green-Naghdi system (3.2) 



'CLi ''"''" "'"'• ^ — ^CL 
up to a remainder R bounded by 

with C = C{M*j^Q,5^,i^i^,5maxi i^maxjin, |A|, \9\), and the corrector term U^ is estimated as 
ll^1L~([o,T*];ff=) + ll^*^1L~([o,T*];ff=) ^ ^ max(e(52 _ -y),^) min(t, Vt). 



Additionally, if there exists a > 1/2, Af^^g, T' > such that 
6 5 

fc=0 i:=0 

then the remainder term, R is uniformly bounded for t G [0, T"] ; 

with C = C(A/^'^g, 5~j„,(5ma2;,/i,„aK, JTi, |A|, \9\) and W^ is uniformly estimated as 

ll^iL~([o,T«];fl--) + ll^*^1L°°([o,T«];H=) - ^ niax(e((52 - 7), /i) min(i, 1). 
Proof. The consistency of the Constantin-Lannes decoupled approximation towards the original 



Green-Naghdi system (4.6) (in the shallow water regime) has been pr oved in |19 . Using the bound- 
it is obvious to check that 



3.1 



edness of U^-^ in relevant spaces, and following the proof of Theorem 

the exact same result holds with regards to our Green-Naghdi system (|3.2p, in the Camassa-Holm 



regime. D 

Remark 8.6. Let us note that Proposition \8.4\ ensures that the above result is not empty, but on the 
contrary is valid for long times (provided that v^' > and the initial data sufficiently smooth). More 
precisely, if the initial data (C°, v^) ^ U^ E H^^'^ , s > sq + 1, sq > 1/2, then for any p G VcH o-nd 
X,9 such that v^' > i/q > 0, then there exists Ci,C2 = C{^imax,S^l„,dmax, |A|, \0\,v^^, |t/°L,+7), 
independent ofp, such that for M*_^q > Ci, the decoupled approximate solution satisfies the uniform 



bound of Proposition 8.5: with T* > Cije. 

Moreover, upon additional condition on the decaying in space of the initial data, there exists 
Ci,C2 such that for any -/Vf'_|_g > Ci, one has T' > C2/max(e, /i). 
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The above properties, thanks to the well-posedness and stabihty of our Green-Naghdi system 
proved in this work, are sufficient to fully justify the solutions of the Constantin-Lannes decoupled 
approximation as approximate solutions of the Green-Naghdi model, and therefore as approximate 



solutions of the full Euler system (2.4) 



Proposition 8.7 (Convergence of the decoupled model). Let p e Vch o,s in (3.1) and A,6' e M 
such that (8.3) holds. Let U^ = (C", ■f/;'')^ G H'^^^ , N sufficiently large, satisfy the hypotheses of 
Theorem 5 in '281 as well as (H1),(H2). Denote by U^ = (CiV')^ ^'ic solu tion of the full Euler 
system (2.4) and by U't,^ the decoupled approximation defined in Definition 8.^ with initial data 



[/°. Then there exists C,T > 0, independent o/p G VcH, such that for any < t < T, one has 

|f^CL-f^lL=^([o,t];if^-) < C {so min(t,ii/2)(i + eot) + Bo-'t), 



with Eq = max(e((5^ — j),iJ,), and U = (Ci "^[C? "0])^ where w[C, V"] ** defined as in (3.3) 



Moreover, if the initial data is sufficiently localized in space, as in the second part of Proposi- 
tion's^ then one has 

l|C^&L-C7lL=o([o,*];if=) ^ ^(^0 min(t,l)(l + e^t) + Bo-^t). 



Proof. By Proposition 

U = Ulr + U" satisfie^T 



VplA such that 



we know that for any Uqi^ , there exists V^ — u \^qi^\ 
le Green-Naghdi system (3.2), up to a small remainder R. Control- 
ling the difference \\U — t/GN*'!!, where I/gn'' is the solution of the Green-Naghdi system (3.2) with 



corresponding initial data, is done exactly as in Theorem |3.4[ and we omit the proof. The result 
is then a straightforward consequence of the triangular inequality, as t/P — C/gn*^ is estimated in 
Theorem 3.4 and ||C^'^[C^cl]|| in Proposition 8.5 D 



A Product and commutator estimates in Sobolev spaces 

Let us recall here some product as well as commutator estimates in Sobolev spaces, used throughout 
the present paper. 

Lemma A.l (product estimates). 

Let s > 0, one has\ff,g e H%R)f]L'^(R), one has 



/ 5 |^= < I / 



9 



H= 



f 



H= 



If s > sq > 1/2, one deduces thanks to continuous embedding of Sobolev spaces, 

I /5 Iff. < I / Iff.l 9 Ins- 
More generally, fors>0 and sa > 1/2, one has V/ G H''{R)C\H''°{R),g G H'{R), 

/I <C I /* I I I _l_ / I /" I I I \ 

9 ifjs r^ \ J \hso I 9 Iffs + \| / \hs\ 9 \hso/s>so ' 

Let F G C°°{R) such that F(0) =0. If g e iJ" (M) f] i°° (K) with s > 0, one has F{g) G H" 
and 

I F{g) 1^, < C{\g\^^,\F\^^}\ n\ 



\Hs 



Finally, we will use 



f 9 f7 = + l ^ C \ f Lr. + l 5 ff^ + l- 

The first estimates are classical (see [l]25 27|), and the last one follows straightforwardly from 

I |2 I |2 I |2 ||2||2 I |2 I |2 

|/ff|ff=+i ^ \ I 9\hs+ l^\dx{f 9)\hs^ \ f \hs\9\hs^ l^\9 dxf \fjs+ l^\ f d^g\^^. 

We know recall commutator estimate, mainly due to the Kato-Ponce pF', and recently improved 
by Lannes [27] (see Theorems 3 and 6): 
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Lemma A. 2 (commutator estimates). 

For any s>0, and d^f,g G L°°(M) f] ^""^W, one has 

I [A',/]3 L2 < \ d^f \hs-i\ 9 \l^ + \ dxf \l^\ g \hs-i- 
Thanks to continuous embedding of Sobolev spaces, one has for s > sq + 1, sq > 2 i 

I [A',/]5 \l2 ^ I d^f Ih.-i| 9 \„s-i- 
More generally, for any s > and sq > 1/2, dxf,g G if'"'(M) f] -ff''^Hlf^)i o"e has 
I [A^/]g 1^. < I d^f I^.J 5 1^.-1+ (I d^f Ih,-i| 9 \hJs>so+i- 
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